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1.  Introduction 


This  report  addresses  the  aerodynamic  forces  and  moments  on  spacecraft  due  to  their  motion  through 
a  rarefied  atmosphere.  The  impetus  for  this  work  arises  from  a  desire  to  capitalize  upon  the  features 
of  an  existing  Aerospace  developed  tool,  ATRIUM,  which  performs  detailed  thermal  analyses  of 
complex  spacecraft  configurations.  The  code  utilizes  randomly  distributed  ray  tracings  in  conjunction 
with  a  suite  of  geometric  primitives  to  calculate  surface  temperature  gradients  due  to  sun  illumination 
and  accounts  for  both  multiple  surface  reflections  and  shadowing.  In  mathematical  terms,  ATRIUM 
provides  a  method  to  perform  surface  integrals  of  incident  rays  over  complex  geometries  and  thus  can 
be  exploited  to  calculate  aerodynamic  forces  and  torques  for  free-molecular  flow. 

On  typical  satellite  configurations  at  altitudes  below  370  km,  the  aerodynamic  torque  frequently 
dominates  [if  and  can  have  a  value  in  the  neighborhood  of  0.0136  Nm.  Although  the  magnitude  of 
this  torque  is  small  by  terrestrial  standards  (as  are  the  other  environmental  disturbances),  it  has  a 
nonnegligible  affect  on  spacecraft  motion  over  an  extended  period  of  time.  Accurate  predictions  of 
the  aerodynamic  disturbances  requires  knowledge  of  the  atmospheric  density,  surface  relative 
molecular  velocity  and  surface  reflection  properties.  Even  at  fixed  altitude  there  will  be  variations  in 
the  density  due  to  Earth’s  diurnal  cycle  and  solar  activity.  To  correctly  compute  the  surface  relative 
molecular  velocity,  it  is  not  sufficient  to  know  the  inertial  motion  of  the  spacecraft;  the  motion  of  the 
atmosphere  is  also  required.  Surface  reflection  properties  are  not  well  known  for  many  materials  and 
experimental  data  indicate  that  they  vary  with  surface  temperature  and  contamination  as  well  as  angle 
of  attack  of  the  incoming  flow  relative  to  the  surface.  In  this  study,  we  adopt  a  simple  gas-surface 
interaction  model  parameterized  by  two  empirical  coefficients  [2].  There  is  a  substantial  literature  on 
this  specific  topic  and  a  variety  of  models  which  have  been  tested  against  experimental  data.  It  should 
be  noted  that  most  density  models  of  the  upper  atmosphere  were  developed  by  observing  the  orbital 
decay  of  satellites  and  are  based  upon  an  assumed  value  of  the  drag  coefficient.  The  drag  coefficient 
however  is  critically  dependent  on  the  nature  of  the  gas-surface  interaction.  It  is  seen  therefore,  that 
proper  modeling  of  the  gas-surface  interaction  is  critical  in  obtaining  accurate  estimates  of  the 
aerodynamic  disturbances. 


At  altitudes  above  120-150  km,  the  atmospheric  density  is  sufficiently  low  that  conventional 
continuum  assumptions  are  no  longer  valid.  A  dimensionless  parameter,  the  Knudsen  number,  is  used 
to  classify  the  flow  regime  along  the  spectrum  from  continuum  to  highly  rarefied.  The  Knudsen 
number,  Kn,  is  defined  by  K n=X/L  where  X  is  the  average  mean  free  path  of  the  molecules  (the 
average  distance  traveled  by  a  molecule  before  collision  with  another  molecule)  and  L  is  a 
characteristic  length  representative  of  the  spacecraft.  By  convention,  a  flow  with  Kn  <  0.01  is 
regarded  to  be  in  the  continuum  regime  where  intermolecular  collisions  dominate.  Flows  with  Kn  > 
10  are  termed  free-molecular  where  molecule-surface  interactions  govern  the  flow.  Intermediate 
values  of  Kn  correspond  to  the  transition  regime  where  molecule-molecule  and  molecule-surface 
interactions  are  of  equal  importance.  Even  for  extremely  low  earth  orbits,  X  is  almost  1  km.  and  hence 
the  calculation  of  aerodynamic  forces  and  torques  upon  the  spacecraft  must  be  based  on  free- 
molecular  flow  theory.  It  is  the  transfer  of  momentum  to  the  surface  by  atmospheric  molecules  which 
determines  the  net  force  and  torque.  Momentum  transfer  occurs  both  when  molecules  arrive  at  the 
surface  and  upon  subsequent  reflection.  The  free-molecular  assumption  allows  us  to  treat  these  cases 
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separately  and  add  their  effects.  Momentum  exchange  due  to  molecules  colliding  with  each  other  are 
neglected  and  only  molecule-surface  interactions  are  important.  In  addition,  geometrically 
complicated  surfaces  can  be  easily  treated  by  partitioning  into  simpler  constituent  surfaces  except  for 
complications  due  to  shadowing  and  multiple  surface  collisions. 
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2.  Nonrotating  Bodies  in  Hyperthermal  Free-molecular  Flow 


As  mentioned  in  the  Introduction,  calculation  of  the  aerodynamic  force  and  moment  on  a  spacecraft 
requires  knowledge  of  the  velocity  of  the  impinging  atmospheric  molecules  relative  to  the  surface.  In 
general,  the  spacecraft  will  be  translating  and  rotating  relative  to  the  atmosphere  and  hence  the 
surface  relative  velocity  of  the  approaching  molecules  will  vary  over  the  surface.  In  addition,  even  for 
a  purely  translational  motion,  the  molecules  will  have  a  distribution  of  velocities  about  a  mean  value 
dependent  on  the  atmospheric  temperature.  In  the  present  section  we  will  assume  that  the  spacecraft 
does  not  rotate  relative  to  the  atmosphere  and  we  shall  also  neglect  the  random  thermal  motion  of  the 
atmospheric  molecules.  The  later  assumption  is  frequently  refereed  to  as  the  hyperthermal 
approximation  and  is  valid  at  sufficiently  low  temperatures  when  the  random  thermal  motion  of  the 
molecules  is  negligible  compared  to  the  speed  of  the  spacecraft.  These  simplifying  assumptions  allow 
us  to  treat  the  incoming  molecules  as  a  collimated  beam  impinging  on  an  element  of  surface  area  with 
a  single  characteristic  velocity  (care  being  taken  to  only  include  those  portions  of  the  surface  exposed 
to  the  oncoming  flow).  The  effects  of  removing  these  assumptions  will  be  addressed  in  subsequent 
sections. 


2.1  Gas-Surface  Interaction 

In  order  to  calculate  the  forces  exerted  upon  the  surface  by  the  impinging  molecules,  we  must  specify 
the  momentum  exchange  that  occurs  between  the  molecules  and  arresting  surface.  By  convention, 
two  canonical  forms  of  surface  reflection  are  defined.  Specular  reflection  (Fig.  2.1)  is  a  deterministic 
process  in  which  the  angle  of  incidence  equals  the  angle  of  reflection  while  the  incident  velocity, 
reflected  velocity  and  surface  normal  are  all  coplanar.  The  speeds  of  the  incident  and  reflected 
molecules  are  not  necessarily  equal  (some  authors  do  however  assume  this  equality  in  their 
definition).  In  the  case  of  diffuse  reflection  (Fig.  2.2),  the  incident  molecules  are  trapped  into  the 
interstices  of  the  surface  losing  all  knowledge  of  incoming  direction  and  subsequently  re-emitted 
from  the  surface  with  a  random  distribution  of  speed  and  random  direction  governed  by  a  “cosine 
distribution”  i.e.  there  are  as  many  molecules  with  a  particular  speed  leaving  the  surface  in  one 
direction  as  leaving  in  the  opposite  direction.  For  both  modes  of  reflection,  the  speed  with  which  the 
molecule  leaves  the  surface  is  determined  by  the  amount  of  energy  it  transfers  to  the  surface  before 
being  re-emitted.  The  degree  of  equilibrium  attained  between  the  molecule  and  surface  is  measured 
by  the  energy  accommodation  coefficient  cx„  defined  by 


Where  Et  ,Er  are  the  incident  and  reflected  energy  fluxes  respectively  and  E„  is  the  energy  flux  that 
would  be  carried  away  by  the  reflected  molecules  if  they  were  in  thermal  equilibrium  with  the 
surface.  Thus  a„=0  corresponds  to  no  energy  exchange  between  the  incident  molecules  and  the 
surface  while  a„=l  characterizes  the  incident  molecules  coming  into  equilibrium  with  the  surface 
before  being  re-emitted. 
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Figure  2.1 .  Specular  Reflection.  Figure  2.2.  Diffuse  Reflection. 


We  next  introduce  the  momentum  accommodation  coefficients  which  relate  the  incident  and  reflected 
momentum  fluxes.  It  is  found  [2]  that  better  agreement  is  attained  between  theoretical  calculations 
and  experimental  data  if  separate  coefficients  are  assigned  for  normal  and  tangential  momentum 
exchange.  Resolving  the  molecular  momentum  into  components  normal  and  tangential  to  the  surface, 
we  define  the  respective  momentum  accommodation  coefficients 


a.  = 


Pi-Pr 
P>~P„ 
r. 


(2.1) 


where  p„pr  are  the  incident  and  reflected  normal  momentum  flux  and  x„Tr  are  the  incident  and 
reflected  tangential  momentum  flux.  The  parameters  pw, xw  denote  the  normal  and  tangential 
components  of  momentum  flux  which  would  be  carried  away  from  the  surface  by  the  diffusely 
reflected  molecules  if  they  were  in  thermal  equilibrium  with  the  surface.  Note  that  tb,=0  by  the 
definition  of  diffuse  reflection.  For  completely  diffuse  reflection,  ct,=1  regardless  of  the  degree  of 
thermal  accommodation  while  o„=l  additionally  requires  complete  thermal  accommodation.  For 
complete  specular  reflection  with  no  energy  exchange  between  the  molecules  and  the  surface, 
CT,=a„=0. 

In  order  to  calculate/?*,,  we  must  find  an  expression  for  the  average  normal  velocity  of  diffusely 
reflected  molecules  which  are  in  thermal  equilibrium  with  the  surface.  With  reference  to  Fig.  2.2,  the 
molecules  have  become  fully  accommodated  to  the  surface  and  leave  with  a  random  speed  v 
characteristic  of  the  surface  temperature  and  at  random  direction  0  (with  respect  to  the  surface 
normal)  governed  by  a  cosine  distribution.  The  kinetic  theory  of  gases  (see  Appendix  A)  provides  us 
with  the  celebrated  Maxwell-Boltzmann  speed  distribution  for  the  molecules  in  equilibrium  with  the 
surface.  If  we  denote  this  speed  distribution  by  0(v),  then  the  probability  that  a  molecule  has  speed  in 
the  range  (v,v  +  dv)  and  direction  between  6  and  &+dd  is  given  by  <p(v)cosddvdd .  The  normal 
component  of  velocity  is  vcosS,  and  hence  the  average  normal  velocity  ,  Vw,  is  given  by 
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nil  °° 

Vw  =  (vcosS)  =  J  j*  v0(v)cos2  QdvdO 

0  0 

Performing  the  integration  over  0  and  expressing  the  mean  speed  in  terms  of  the  surface  temperature 
TW[A.17,  A.18]  we  obtain  Vw  =  (nkTw/2m)'12  where  k  is  the  Boltzmann  constant  (1.38065  xl0~23 
joule/molecule  °  K)  and  m  is  the  mass  of  a  single  gas  molecule.  Recalling  that  k  =  "R/N0  where  %  is 
the  universal  gas  constant  (8314.5  joules/kg-mole°K)  and  N0  is  Avogadro's  number  (6.02214  xlO26 
molecules/kg-mole),  we  can  express  Vw  in  the  form 


V  = 


(2.2) 


where  7K  denotes  the  molecular  weight.  It  sometimes  proves  convenient  to  express  the  above  equation 
in  terms  of  the  gas  constant  R,  which  unlike  *R,  depends  upon  the  particular  gas  of  interest.  In  terms  of 
the  universal  gas  constant,  the  equation  of  state  of  an  ideal  gas  is  pV  =  NKT / N0  where  p  is  the 
pressure,  V  is  the  volume,  and  N  denotes  the  total  number  of  molecules.  Comparing  this  with  the 
alternative  form  p  =  pRT  where  p  is  the  gas  density,  we  find  that  R-KM.  Hence  eq.(2.2)  can  be 
expressed  in  the  alternate  form 


(2.3) 


2.2  Aerodynamic  Forces  and  Moments 

Our  first  model  is  developed  under  the  assumptions  that  1)  the  body  is  not  spinning  relative  to  the 
atmosphere  and  2)  the  mean  thermal  motion  of  the  atmosphere  is  negligible  compared  to  the  speed  of 
the  body  through  the  atmosphere;  the  so  called  hyperthermal  flow  assumption. 

Figure  2.3  shows  atmospheric  molecules  impinging  upon  the  spacecraft  surface  5  which  is  assumed 
to  be  convex  so  as  to  avoid  complications  due  to  shadowing.  The  molecules  all  have  velocity  v 
relative  to  the  surface  so  that  the  incident  flow  on  an  element  of  surface  area  dA  can  be  viewed  as  a 
collimated  molecular  beam  of  mass  density  p.  Denote  by  nin  a  unit  inward  drawn  normal  to  the 
surface  and  let  a  be  the  angle  between  hm  and  V. 

cosa  =  v-nia  (2.4) 

Where  v  =  Vy/|v|.  It  is  important  to  note  that  no  impingement  occurs  on  those  portions  of  the  surface 

where  cos  a  <  0 .  Hence  we  will  insert  the  factor  H(cosa)  (H  denotes  the  Heaviside  step  function) 
into  the  surface  integrals  below  to  effectively  limit  the  aerodynamic  pressure  to  those  portions  of  the 
surface  which  actually  experience  impingement.  For  an  arbitrary  surface  S  as  shown  in  Fig.  2.3,  there 
is  no  natural  set  of  basis  vectors  to  resolve  the  shear  stress  on  the  element  of  area  dA  so  we  will 
establish  these  vectors  with  respect  to  the  direction  v  of  the  incident  molecules. 


5 


\ 


Figure  2.3.  Spacecraft  Surface  in  Free-Molecular  Flow. 


Decomposing  v  into  a  component  normal  to  the  surface  and  its  orthogonal  projection  into  the  tangent 
plane,  we  can  write 

v  =  cosaiiijn  +  sin  at  (2.5) 

Which  serves  to  define  the  unit  tangent  direction  t.  Forming  the  triple  cross  productnin  x(vxnin), 
yields  the  following  expression  for  t  directly  in  terms  of  the  unit  vectors  V  and  njn 

t_  fiinx(vxn.n) 

II **j„  x(vx nin ) || 

The  vectors  |t,iiin,  t  x  njn|  constitute  an  orthonormal  triad  on  the  surface. 

In  a  time  interval  At  the  mass  delivered  to  dA  is  contained  in  a  parallelepiped  of  base  area  dA  and 
slant  height  VA  t  which  constitutes  a  volume  VAf  cos  a  dA.  Thus  the  mass  flux  SQ  is  given  by 


5Q  =  pV  cos  a  (2.6) 

With  the  aid  of  eq.(2.5),  the  incident  momentum  flux  V  SQ  resolves  into  the  components 


p.  =  V  cosaSQ 
t,  =  Vsma5Q 


(2.7) 


Under  steady  state  conditions  the  incident  and  reflected  mass  fluxes  are  equal;  hence  pw=Vw  SQ 
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where  Vw  is  given  by  eq.(2.3).  The  reflected  momentum  flux  is  given  by  the  expression:  — pr Sin  +  Tr  t. 

Subtracting  the  reflected  from  the  incident  momentum  flux  yields  the  force  per  unit  area,  dt/dA, 
exerted  by  the  molecules  upon  the  surface. 

^  =  {pi  +  Pr)nm  +  (T,  -  *,)t  (2.8) 

The  two  components  of  reflected  momentum  flux  can  be  expressed  in  terms  of  the  incident  fluxes,  p„, 
and  momentum  accommodation  coefficients  through  eq.(2.1)  By  utilizing  eqs.(2.6)  &  (2.7),  we  can 
write  eq.(2.8)  in  the  form 

d(  ,  IT,  s  V  1-  J 

—  =  pV  cosaj  (2  -  an  -  a, ) cos  a  +  cx„y  ni„  +  j-  (2.9) 

where  we  have  eliminated  t  in  terms  of  niuand  v  through  eq.(2.5).  This  form  is  simpler  to  integrate 

since  v  (unlike  t)  will  not  vary  over  the  surface.  Integrating  over  those  portions  of  the  surface  which 
are  impinged  by  the  atmospheric  molecules,  and  assuming  that  the  atmospheric  density, 
accommodation  coefficients,  and  surface  temperature  are  constant  over  the  surface,  we  obtain  the  net 
force  f  and  moment  g  about  point  O  (see  Fig. 2.3).  All  vectors  are  assumed  to  be  resolved  in  a  set  of 
body-fixed  axes  and  indicated  by  bold  face  type. 

(  =  PV2[o,Atv  +  oXvjV)T,+(2-a„ -a,)n]  (2.10) 

%  =  pV![ct,A,c,  xit  ajvjvjg  +  (2-ct, -cr,)g]  (2.11) 

where  the  newly  introduced  terms  are  only  dependent  upon  the  surface  geometry  and  angle  between 
the  incident  flow  and  surface  normal. 

AP  =  H(cos  a)  cos  a  dA  Apcp  =^//(cosa)cosardA 

7f=|j>//(cosa)cosanta  dA  !7^  =  ££//(cosa:)cos2a  nin  dA  (2.12) 

^  =  ^ //(cos  a)  cos  a  rxnhi  dA  =^H(cosa)cos2a  rxn  indA 
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2.3  Closed  Form  Solutions 

As  a  means  of  checking  the  implementation  of  the  aerodynamic  force  and  torque  calculation,  we 
provide  closed-form  solutions  for  the  following  geometric  shapes. 

•  Sphere 

•  Flat  Plate 

•  Right  Circular  Cylinder 

•  Circular  Cone  Frustum 

The  force  and  moment  components  are  plotted  as  functions  of  the  angle  of  attack  (except  for  the 
sphere) 

2.3.1  Sphere  in  Hyperthermal  Free-molecular  Flow 

Figure  2.4  shows  a  sphere  of  radius  a  with  a  system  of  spherical  coordinates  { r,6,(p }  at  its  center. 
Without  loss  of  generality,  we  may  assume  that  the  incident  molecular  flow  V  is  directed  along  the  z 
axis.  The  surface  of  the  sphere  is  described  by:  r=a ,  0 <Q<n,  0<(p<2n  with  element  of  surface 
area  dA  =  a 2  sin  6dQd(p.  The  unit  inward  normal  is  given  by 

Hi,,  =  -sin0cos^>i-sin0sin<p  j-cos0k 

where  ij,k  denote  unit  vectors  along  the  axes  of  x,y,z  respectively.  Recognizing  that  the  factor 
H(cosa )  limits  the  integration  to  n/2<B<n,  the  surface  integrals  are  readily  computed  and  have 
the  values:  Ap=7ta2,  Tl=2na2l3  k,  ?r2=7W2/2  k  .  Hence  the  net  force  on  the  sphere  is  a  drag  with 
magnitude  D  given  by 

D  =  ±na2pV2{2  +  o,-an+^onty  (2.13) 

It  is  clear  from  symmetry  (and  verified  by  calculation)  that  the  aerodynamic  torque  about  the  center 
of  the  sphere  is  zero. 


Figure  2.4.  Sphere  in  Free-Molecular  Flow. 


8 


2.3.2  Flat  Plate  in  Hyperthermal  Free-molecular  Flow 

Figure  2.5  shows  a  flat  plate  of  area  A  at  angle  of  attack  P  with  respect  to  the  incident  molecular 
velocity  V.  A  rectangular  coordinate  system  is  established  with  origin  at  the  centroid  of  the  plate  with 
x,y  axes  respectively  parallel  and  perpendicular  to  V.  Here  cos  a  =  sin  /3  >  0  for  angles  of  attack 
between  zero  and  180  degrees.  Since  the  surface  is  flat,  the  integrals  are  trivial  to  evaluate. 


Ap=A  sin/3,  9{  =  Asinf}nm,  %  =  A  sin2/3  nin 

The  resolution  of  the  inner  surface  normal  (with  respect  to  the  direction  of  incident  flow)  along  the 
drag  ( x )  and  lift  (y)  directions  is  given  by  nbi  =  sin /3i  +  cos/3  j.  The  net  force  on  the  plate  then 
assumes  the  form 

f  =  ipV2A(CDi  +  Cj) 


Where  the  lift  and  drag  coefficients  are  given  by 


Ct  = 


^  -^7  +  (2  ~  -crr)sm/3 


sin2j6 


Cd=2 


<*,  +  an  y sin  P  +  (2  -  -  a < )  sin2)3 


sin/3 


(2.14) 

(2.15) 


It  is  readily  shown  that  the  aerodynamic  moment  about  the  centroid  of  the  plate  is  zero. 

The  lift  and  drag  coefficients  are  plotted  below  as  functions  of  the  angle  of  attack  for  cr„=  a, =0.7  and 
Vw/V= 0.05.  In  general,  the  drag  attains  its  maximum  value  at  P=7t/2  being  symmetric  about  this  point 
and  vanishing  when  P=0,n  while  the  lift  is  anti  symmetrical  about  P=ti/2  vanishing  at  the  points 
P=0,7t/2,7t. 
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2.3.3  Right  Circular  Cylinder  in  Hyperthermal  Free-molecular  Flow 


► 

x 


Figure  2.7.  Cylinder  in  Free-molecular  Flow 


Figure  2.7  shows  a  right  circular  cylinder  of  length  l  and  base  radius  a  at  angle  of  attack  P  with 
respect  to  the  incident  molecular  velocity  V.  A  rectangular  coordinate  system  with  origin  at  the 
centroid  O  of  the  cylinder  has  its  2  axis  along  the  axis  of  the  cylinder  and  x  axis  in  the  plane 
determined  by  Y  and  the  cylinder  axis.  It  suffices  to  consider  angles  of  attack  in  the  range 
0<  P<7i  /  2.  For  purposes  of  computing  the  surface  integrals,  we  partition  the  cylindrical  surface 
into  its  end  bases  S,,S2  and  lateral  surface  S3. 

•  Force  &  Moment  on  Lateral  Surface 

In  terms  of  the  cylindrical  coordinates  (r,6,z)(x  =  rcosd,  y  =  rsin0),  the  surface  S3  is  described  by 
r  =  a,  0  <  6  <  2n,  —If 2  <  z  <  Z/2  and  the  element  of  surface  area  is  given  by  clA  =  addclz  ■  Observing 
that  iijj,  =  -cos0i  -sin0  j  and  v  =  sin/?i  +  cos/3k.  we  have  cosa  =  -sin/3cos0;  thus  the  region  of 
integration  will  be  limited  to  n  /  2<6  <3n  /  2  for  the  specified  range  of  angle  of  attack. 
Straightforward  integration  yields  the  following  expressions  over  the  lateral  surface. 

Afl  =  2alsmp  Apcl)=-7t/2a2lsml3i  ^  =  0 
7f=7r/2n/sin/3i  %  =  4/3a/sin2/I  i  =  ® 


11 


Resolving  the  force  and  moment  along  the  {x,y,z}  axes,  we  find  the  force  f3  and  moment  g3  (about 
point  O)  on  the  lateral  surface. 


f3  =  pV-al< 


n 


V  (8  4  2  \  .  0 

cr  —  + - <7+  —  cr,  sin  >3 

V  V3  3  3  J 


sin  P  i  +  tr,  sin  2)3  k 


g3  =  tt/4  pV/2rt2/crf  sin  2>S  j 


(2.16) 

(2-17) 


•  Force  &  Moment  on  End  Surfaces 

The  end  face  S2  of  the  cylinder  is  not  impinged  upon  by  the  flow  over  the  specified  range  of  angle  of 
attack.  To  calculate  the  force  and  torque  on  5,,  we  apply  the  results  obtained  above  for  the  flat  plate 
with  A- tui2  (see  eqs.(2.14)  &  (2.15))  Noting  from  Fig. 2.7  that  the  angle  of  attack  of  5,  with  respect  to 
V  is  P+rc/2,  and  resolving  the  drag  and  lift  forces  along  the  coordinate  directions  {a-,v.z},  we  find  that 
the  force  f,  on  5,  is  given  by 


f,  =7r/2aJpV,3jcr(sin2j3i  +  2cos/3 


+  (2-a„)cosP 


(2.18) 


The  aerodynamic  moment  on  *S,  about  point  O  is:  -/  /  2k  x  f,  which  can  be  shown  to  be  equal  to  -g3. 


•  Net  Force  &  Moment  on  Cylinder 

Adding  the  contributions  from  »S|  and  we  find  that  the  net  force / on  the  cylinder  can  be  written  in 
the  form  f=fx  i  +/.  k  where 


f^jalpV1  sin  13 
6 


/  =  alpV 2  cos  (i 


4(4  +  o,  -  2a„ )sin  /?  +  6n  j  o,  cos  p  +  3nan 
2a,smP  +  n^(2-a„)cosp  +  n^an 


(2.19) 


The  net  moment  about  the  centroid  of  the  cylinder  is  zero  at  any  angle  of  attack. 

The  two  force  components  are  plotted  below  as  functions  of  the  angle  of  attack  for  cr,  =  cr,  =  0.7 
Vw/V= 0.05,  and  a/l= 0.2. 
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Forces  on  Cylinder 


2.3.4  Cone  Frustum  in  Hyperthermal  Free-molecular  Flow 

Figure  2.9  shows  a  cone  frustum  of  circular  cross  section  in  a  free-molecular  flow  where  the  incident 
velocity  V  is  at  an  angle  of  attack  P  with  respect  to  its  axis  of  symmetry.  The  end  surfaces  S,,S2  of  the 

body  are  circular  disks  of  radii  RhR2  respectively  and  y  denotes  the  half  angle  of  the  cone.  We  restrict 
our  attention  to  the  range  of  angles  0  <  y  <  n/2,  0  <  P  <  tc/2  .  A  rectangular  coordinate  system  {x,y,z} 
is  introduced  with  origin  at  the  apex  of  the  cone  and  z  axis  coincident  with  the  axis  of  symmetry. 
Without  loss  of  generality,  we  choose  the  x  axis  such  that  V  lies  in  the  x-z  plane.  In  terms  of 
cylindrical  coordinates  {r,0,z},  the  lateral  surface  S3  of  the  cone  frustum  is  described  by:  0<6  <2n , 
z,  <  z  <  z2  where  z;  =  Rt  coty  ( i  =  1,2)  and  r{z)  =  tan  y  z-  On  S3,  the  position  vector  to  an  arbitrary 
point  is  given  by  r-  tanyzcos0i  +  tanyzsinSj  +zk  and  the  element  of  surface  area  is 
dA  =  siny/cos2  y  zdzdO.  From  the  figure,  nin  =  -cosy  cos0  i  -  cosy  sin0  j  +  siny  k  on  S3  and 
v  =  sin  (3  i  +  cos  [i  k .  Forming  the  inner  product  of  these  two  vectors  we  find 

cos  a  =  cos  P  sin  y  -  sin  P  cos  y  cos  6  (2.20) 

To  determine  the  portions  of  S3  that  are  impinged  by  the  atmospheric  molecules,  we  need  to 
determine  the  sign  of  cosa  for  0<@<27r.  Since  cosa  vanishes  when  cos0  =  tany/tan P,  we 
consider  separately  the  two  cases:  1)  0  <  P  <  y  and  2)  y  <  P  <  n/2. 
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Figure  2.9.  Cone  frustum  in  Free-molecular  Flow. 


•  Case  1 :  0  <  P  <  y 

It  is  clear  from  eq.(2.20)  that  in  this  case  cosa  >0  for  all  values  of  6,  hence  the  entire  lateral  surface 
of  the  cone  intercepts  the  molecular  flux.  Performing  the  integrations  indicated  in  eq.(2.12)  we  obtain 

Ap=Jt(/$-R*)casP 

7\=^-(R22-  R2  )(cot y  cos y  sin  P i  +  2  sin y  cos P  k) 

7C  1 

%  =  —  (/?22  -  R2)  cos2y  sin2/3i  +  ~(3  +  cos2/?-cos2y  -3cos2/J  cos2y)k 
Apcv  =-^coty(y?j  -  Z?,3 )(—  sin  /3  i  +  2  cos  /3  k) 

^ [R]  ~  R2 ) sin  J3  cot  y  esc  y  j 
&  =  j(/^-/?,,)sin2iScoty  j 

The  force  f3  and  moment  g3  (about  point  O)  on  the  lateral  surface  follow  from  eqs.  (2.10)  &  (2.11). 

f ,=/  o)  sin  2/3  +  crn cot /cosy  sin /?  + (2 -<7„ -<7,)cos2y  sin2/3  i  + 

f  [la,  cos2 p  +  2a„-^sinycos/3  +  (2.21) 

(2  -  crn  -  a,  )(3  +  cos  ip  -  cos  2y  -3  cos  2y  cos  2)3)]  k 
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(2.22) 


g 3  =jpV2(/?2  ~/?f)  cr„-^-cscy +(4-2o-(i +(t,)cosj3  cotysin/3j 

where  /  =  —  pV27i(/?22  -  i?,2) 

We  note  that  in  the  limit  as  y  — »tz:/ 2,  the  lateral  surface  degenerates  into  an  annular  disk  and  the 
above  expressions  are  in  agreement  with  those  obtained  previously  for  a  flat  plate  (see 
eqs.(2.14),(2.15))  corresponding  to  the  area  A  =  ti[r\  -  Rf ). 

•  Case  2:  y  <  (3  <  nj 2 

To  establish  the  0  limits  of  integration  in  this  case,  we  have  plotted  in  Figure  2.10  the  function  cosa 
of  eq.(2.20)  for  6  over  the  interval  (0,2tt).  We  note  that  this  function  is  negative  at  0=0,  positive  at 
0=7T ,  and  symmetric  about  the  latter  point.  The  root  O'  of  the  function  in  the  interval  (0,7t)  is  given  by 

8 *  =  cos-1  {tan  y /tan  /?)  (2.23) 

hence  the  integration  limits  in  this  case  are  0*  < 0 < 2n -0*; z]  <z^z2- 


Figure  2.10.  Plot  of  cos  a(d)  for  (3>y  ■ 


We  will  have  a  continuous  transition  from  Case  1  since  0*  — >0as  /J— »y+.  Performing  the 
integrations  indicated  in  eq.(2.12)  we  obtain 

Ap  =(/?22  -/?2)[(7r-0*)cos/3  +  cotysin0’sin^] 
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Combining  the  results  from  the  above  two  cases  allows  us  to  compute  the  aerodynamic  force  and 
moment  on  the  lateral  surface  of  the  cone  frustum  over  the  full  range  of  angle  of  attack  (0  <  (3  <  njl') 
as  illustrated  in  the  two  plots  below.  The  numerical  results  were  generated  for  a  cone  half-angle  of  60 
deg..  a,=oy=0.7,  and  Vw/V= 0.05.  The  force  and  moment  were  rendered  dimensionless  by  dividing  by 
the  respective  scale  factors  1/2  pV2n^  -  and  1/3 pV2 .  To  find  the  net  force  and 

moment  on  the  cone  frustum,  we  simply  add  the  contribution  from  the  end  surface  S,  which  is  a  flat 
plate  of  area  n Rf  (see  eqs.  (2.14)  &  (2.15)). 
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Figure  2.11.  Forces  on  Cone  Frustum  in  Hyperthermal  Free-molecular  Flow. 


Figure  2.12.  Torque  on  Cone  Frustum  in  Hyperthermal  Free-molecular  Flow. 
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3.  Nonrotating  Bodies  in  Free-molecular  Flow 

In  the  previous  section,  the  model  for  predicting  the  aerodynamic  disturbances  on  a  body  in  a  free- 
molecular  flow  was  based  upon  the  assumption  that  the  body  was  undergoing  a  pure  translation 
relative  to  the  atmosphere  and  that  the  random  thermal  motion  of  the  atmospheric  gas  molecules  was 
negligible  in  comparison  with  the  speed  of  the  body  relative  to  the  atmosphere.  In  this  section  we 
develop  expressions  for  the  aerodynamic  force  and  torque  which  include  the  effects  of  the  random 
thermal  motion  of  the  gas  molecules  (the  body  is  still  assumed  not  to  rotate  relative  to  the  mean 
motion  of  the  atmosphere).  As  will  be  seen,  the  quadratures  which  are  called  for  in  calculating  the  net 
force  an  torque  are  much  more  complicated  that  those  derived  previously  under  the  hyperthermal 
flow  assumption. 


3.1  Aerodynamics  of  a  Body  Translating  in  a  Maxwellian  Flow 

If  we  consider  a  gas  in  equilibrium,  then  the  number  of  molecules  with  random  velocities  between  u 
and  u+du  is  given  by  [A.  19] 


N 


n 


3/2 


mN0 

2“RT 


,3/2 


exp 


mN, 


2KT 


—  n2 


dui  du2  du3 


where  N  is  the  total  number  of  molecules,  m  is  the  mass  of  a  single  gas  molecule  and  T„  is  the 
atmospheric  temperature.  The  constants  *R  and  N0  denote  the  universal  gas  constant  and  Avagadro's 
number  respectively.  Dividing  this  expression  by  the  gas  volume  and  recognizing  the  product  mN0  as 
the  gas  molecular  weight  7%,  we  find  that  the  number  of  molecules  per  unit  volume  with  random 
velocities  between  u  and  u+du  can  be  expressed  in  the  form 


P 

mvyn 


exp^ 


dul  du2  du3 


where  p  is  the  gas  density  and  the  parameter  Va  represents  the  most  probable  speed  [A.16]. 


(3.1) 


V  = 


(3.2) 


The  random  velocity  u  appearing  above  is  with  respect  to  the  mean  motion  of  the  gas.  Denoting  by  V 
the  mean  velocity  of  the  gas,  and  resolving  along  the  unit  triad  {nin,t,njn  xt}  at  the  element  of 
surface  area  dA  (see  Fig.  3.1),  we  can  write 


V  =  V(cosanjn  +  sin  at). 

It  follows  that  the  velocity  v  of  the  molecules  with  respect  to  the  surface  is  given  by  v  =  V  +  u. 
Writing  v  =  v,  +  v2t  +  v3  nin  x  t ,  we  can  express  the  square  of  the  random  speed  in  the  form 

u 2  =(v-V)-(v-V)  =  v2  +V2  -2V(v,cosa  +  v2 sina) 


Hence  the  number  of  molecules  per  unit  volume  whose  surface  relative  velocity  lie  between  v  and 
v+dv  is  given  byy(v,,v2,v3)  dV]dv2  dv3  where 
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V 


Figure  3.1.  Unit  Triad  on  Surface 


/(VI»V2V3) 


pexp(-S2) 

mY>w 


'IS,  .  x‘ 

exp 

1  K,2J 

exp 

—  (v,  cosa  +  v2  sin  a) 

(3.3) 


and  S  denotes  the  molecular  speed  ratio 


(3.4) 


Our  first  model,  based  upon  the  hyperthermal  flow  assumption,  corresponds  to  S  =  <*>  (in  the  sense 
that  V=0). 


Consider  a  stream  of  molecules  with  velocities  in  the  range  (v,v+dv)  incident  upon  the  element  of 
surface  area  dA.  In  a  short  time  interval  At,  these  molecules  sweep  out  a  volume  dAv^At.  Thus 
f(\)dyldv2dvi  dA  v,  Atm  is  the  mass  delivered  to  dA  in  time  At  from  molecules  with  surface  relative 
velocity  in  the  range  (v,v+dv).  The  net  mass  flux  8Q  from  all  incident  molecules  is  given  by 


5Q  =  mj  J  J/(v1,v2,v,)v1rfv1rfv2£?v3 
0 

Note  the  integration  limits  on  v,;  a  molecule  with  v,  <  0  is  traveling  away  from  the  surface.  Utilizing 
eq.(3.3)  and  performing  the  integration,  we  obtain 


where 


52  =  pVari(5cos«) 

f,  (x)  =  ~j=  [exp(-x2 )  +  sfnx(l  +  erf  x)] 


(3.5) 

(3.6) 


The  incident  momentum  flux  due  to  molecules  with  surface  relative  velocities  in  the  range  (v,v+dv), 
is  given  by 
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'n/Mv^v.n.,  +  v2t  +  v,nill  xt )dv]dw2d\i 

Integrating  over  the  region  {0<v,  <oo -oo<v2  <  oo,-oo  <  v3  <«=}  and  noting  that  /(v)  is  an  even 
function  of  v3  (so  that  the  last  component  integrates  to  zero)  we  arrive  at  the  following  expressions  for 
the  components  of  the  incident  momentum  flux  along  the  normal  and  tangential  directions  to  the 
element  of  surface  area. 

A=pV/2r2(5cos«) 

T,  =  pV"S  sin  a  T,  (5  cos  a) 

where  r2  (x)  =  [*  exp(-x2 )  +  /2  j(l  +  2x2  )(l  +  erf  x)  j 

The  functions  F,(x)  and  r2(x)are  plotted  below  in  Figure  3.2 


(3.7) 

(3.8) 


Under  steady  state  conditions',  the  incident  and  reflected  mass  fluxes  are  equal;  hence  pw=Vw  5Q 
where  Vw  and  8<2  are  given  by  eqs.(2.2)  &  (3.5)  respectively.  The  force  per  unit  area  is  given  by 
eq.(2.8),  where  again,  the  two  components  of  reflected  momentum  flux  can  be  expressed  in  terms  of 
the  incident  fluxes,  pw,  and  momentum  accommodation  coefficients  through  eq.(2.1).  Utilizing 
eqs.(3.5)  &  (3.7)  we  finally  arrive  at  the  result 


*L 

dA 


=pv: 


(2  -  cr„  )r2  (5  cos  a)  -  <t,S  cos  a  F,  (5  cos  a)  +  an  -  T,  (5  cos  a) 


n,„  + 


(3.9) 


pV] <7,5F,(5  cos  a)  v 


where  we  have  eliminated  t  in  terms  of  nin  and  V  through  eq.(2.5)  (recall  that  v  =  V/||V||) 


Figure  3.2.  Universal  Functions  T,  and  r2- 
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In  contrast  to  the  hyperthermal  case  treated  earlier,  all  portions  of  the  (convex)  surface  experience 
impingement  due  to  the  random  motion  of  the  molecules;  hence  we  do  not  restrict  the  integration  to 
those  portions  of  the  surface  where  cos  a  >  0 .  As  would  be  expected,  the  molecular  flux  is  maximum 
at  a=0  and  minimum  at  a=n  (see  eq.(3.5)  and  Fig.  3.2).  Integrating  eq.(3.9)  over  the  entire  surface 
and  assuming  that  the  atmospheric  density,  accommodation  coefficients,  and  surface  temperature  are 
constant,  we  obtain  the  net  force  /in  the  form 


In  performing  the  surface  integrals,  recall  that  a  is  the  angle  between  the  surface  unit  inner  normal  nin 
and  mean  flow  velocity  V. 

As  remarked  above,  the  hyperthermal  flow  assumption  (which  neglects  the  random  motion  of  the 
atmospheric  molecules)  corresponds  to  an  infinite  value  of  the  molecular  speed  ratio  S.  Recalling  the 
limiting  values  of  the  error  function:  lim  erf(x)  =  ±1,  it  is  easy  to  show  that 

A— *±“> 


Hm[i  +  erf(5cosa)j  =  2 //(cosa ) 

..  r(Scosa)  . 

hm— - -  =  cosa  H(  cosa) 

5 


From  eqs.(3.6)  &  (3.8)  it  can  be  shown  that 


and  hence 


r2(x)  =  xr,(x)  +  -(l  +  erfx) 

r,(Scosa)  ,  ...  v 

lim— — - -  =  cos'aH(cosa) 

S 


(3.14) 


(3.15) 


Armed  with  these  results,  it  is  readily  demonstrated  that  in  the  limit  of  infinite  S,  eqs.(3.10)  &  (3.12) 
do  indeed  recover  the  hyperthermal  results  given  by  eqs.(2.10)  &.  (2.11). 


21 


3.2  Closed  Form  Solutions 

Although  the  integrations  are  rather  involved,  we  can  obtain  closed  form  solutions  for  the 
aerodynamic  forces  and  torques  for  finite  molecular  speed  ratio  in  the  case  of  simple  geometries. 
Below  we  present  results  for  the  a)  Sphere,  b)  Flat  Plate  &  c)  Right  Circular  Cylinder.  The  force  and 
moment  components  are  plotted  against  angle  of  attack  (except  for  the  sphere)  for  selected  values  of 
the  molecular  speed  ratio  and  are  compared  with  the  hyperthermal  results.  We  demonstrate 
numerically,  and  prove  analytically,  that  the  resulting  force  and  moment  expressions  reduce  to  the 
previously  obtained  hyperthermal  results  as  the  molecular  speed  ratio  tends  to  infinity. 

3.2.1  Sphere  In  Free-molecular  Flow 

Figure  2.4  shows  a  sphere  of  radius  a  with  a  system  of  spherical  coordinates  {r,6,(p)  at  its  center. 
Without  loss  of  generality,  we  may  assume  that  the  incident  mean  molecular  flow  V  is  parallel  to  the 
z  axis.  The  surface  of  the  sphere  is  described  by:  r=a ,  0  <9<n,  0<(p<2n  with  element  of  surface 
area  dA  =  cr  sin  Bdddcp.  The  unit  inward  normal  is  given  by 

n:„  =  -sinScos  i-sin0sin<p  j-cos0k 

where  i,j,k  denote  unit  vectors  along  the  axes  of  x,y,z  respectively. 


To  evaluate  the  force  associated  surface  integrals  appearing  in  eq.(3.11)  we  make  the  substitution 
x  =  -Scosd .  The  integration  over  <p  is  immediate  and  we  find 

F,  =  Ina1  / S J  ri(jc)<k 
Ta=2m,2/S2  fsxr,(x)dxk 
T,  =  2m2/S'fsx2ri(x)dxk 
%  =  2m2 /S2  I*  jcr2(jt)drk 

J-S 


Utilizing  integration  by  parts,  we  can  perform  the  integration  over  A'. 

S  l  i 

J  r,  (x)dx  =  -  (l  +  252  )erf (5)  +  Se 's'’ 

-5  ^ 

{  x%  (x) dx  =  (l  +  4S4 )  erf (5)  +  —  (252  -  l)5e-5’’ 

J  X r2(x)dx  =  j^(454  +  4 S2  -  l)erf(S)  +  -^~(2S2  +  l)5e's! 

Combining  the  above  results  we  find  that  the  net  force  on  the  sphere  is  a  drag  with  magnitude 
pV2na 2  CDj 2  where  the  drag  coefficient  is  given  by 


(3.16) 


CD  = 


2  -  cr  +  cr. 


2  S'1 


4S4  +4S2  -1  252  + 1 

•erf (S)+ — 7=—  e  s 
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4  V 
+  -  <T„  — 

3  "  V 


(3.17) 
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Taking  the  limit  of  this  expression  as  5— >  we  recover  the  hyperthermal  result  in  eq.(2.13). 

The  drag  coefficient  is  plotted  below  as  a  function  of  S  with  c,=a,= 0.7  and  Vu/V=0.05.  Since  the  drag 
coefficient  decreases  with  S,  we  conclude  that  the  drag  force  is  underestimated  by  the  hyperthermal 
assumption. 


Drag  Coefficinet  on  Sphere 


Figure  3.3  CD  for  Sphere  as  a  Function  of  S. 

It  is  clear  from  symmetry  (and  verified  by  calculation)  that  the  aerodynamic  torque  about  the  center 
of  the  sphere  is  zero. 
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3.2.2  Flat  Plate  in  Free-molecular  Flow 


Figure  2.5  shows  a  flat  plate  of  area  A  at  angle  of  attack  (3  with  respect  to  the  incident  mean  molecular 
velocity  V.  A  rectangular  coordinate  system  is  established  with  origin  at  the  centroid  of  the  plate  with 
x.y  axes  respectively  parallel  and  perpendicular  to  V.  In  contrast  to  the  prior  hyperthermal  analysis, 
the  random  motion  of  the  molecules  will  result  in  impingement  upon  the  back  surface  of  the  plate  and 
hence  we  must  now  also  integrate  over  that  surface.  Denote  by  S+  and  S'  the  front  and  back  surfaces 
(relative  to  V)  of  the  plate  respectively.  The  angle  a  assumes  constant  values  over  the  two  surfaces, 
so  the  surface  integrals  are  trivial  to  evaluate. 

On  S+,  cosa  =  sin/3  and  the  unit  inward  normal  is  given  by  n  =  sin /3i  +  cos/3 j.  Denoting  the 
contribution  to  the  force  (eq.)  (3.10)  by  F.  we  find 


r  _  PV2A 


<T,r,  (Ssin  /3) i  +  ~~JL ( S sin  /3) n  -  <T,r, (S sin  0) sin  /3  n  +  <7„  ~  f, (S sin  /?) n 


On  S  ,  cosa  = -sin/3  and  the  unit  inward  normal  is  (-n).  The  corresponding  force  contribution  F  is 
given  by 


pV2A 


tr,r,  ( -S sin  /3)  i - ■^JL  r2  (-S sin  /3)  n  -  <7,r,  (-5 sin  (3)  sin  p  n  -  crn  —  F,  (-5 sin  0)  n 


The  following  relations  can  be  immediately  verified  and  will  prove  useful  in  simplifying  the 
expression  for  the  net  force  on  the  plate  (see  eqs.(3.6)  &(3.8). 


T,  (x)  +  T,  (-a)  =  — e"A"  +  x  erf  (a) 

■\In 

r,(A)-r,(-A)=A 

r:  (a)  -  T2  (-a)  =  xe~*  +^-(l  +  2A2)erf(A) 


Adding  the  two  force  components  F  and  F  and  introducing  the  expression  for  the  unit  normal  vector 
n,  we  obtain  the  net  force  f  on  the  plate  in  the  form 

f  =  ipFM(CDi  +  Cj) 

where  the  lift  and  drag  coefficients  are  given  by 


Cl=  -2— °±exp(-52  sin2  jg)  +  cr„  sin  2/3  + 

-  O  _ 

— TJL  +  2(2-an  -CTf)sin2/8  cosj3erf(Ssin j3) 

iS 


(3.18) 
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CD  =  2cr„ -^sin2/3  +  -^=- [(2 -crn) sin 2 /3  +  <7,  cos2/?  ]exp(-S2  sin2/3)  + 


(2  -  <r„  f  sin2/3  +  — 1  +  o',  cos2/3 


(3.19) 


sin/lerf(Ssin  jB) 


Taking  the  limit  of  the  above  expression  as  S— we  recover  the  hyperthermal  results  of  eqs.(2.14) 
&  (2.15). 

It  is  readily  shown  that  the  aerodynamic  moment  about  the  centroid  of  the  plate  is  zero. 

The  lift  and  drag  coefficients  are  plotted  below  as  functions  of  the  angle  of  attack  for  <Jn=o,= 0.7  and 
Vw/V=0.05.  The  solid  curves  correspond  to  5=3  while  the  dashed  curves  are  the  corresponding 
hyperthermal  results.  As  in  the  hyperthermal  case,  the  lift  is  antisymmetric  about  P=Jt/2  vanishing  at 
(3=0,rt/2  and  n.  Note  however  that  the  angle  of  attack  corresponding  to  maximum  lift  is  affected  by 
the  molecular  speed  ratio.  The  drag  remains  symmetric  about  the  point  P=rt/2  at  which  it  is 
maximum;  however  it  assumes  a  nonzero  value  at  P=0,7t  for  finite  S  and  (7t  A  0 .  We  deduce  from  the 
plots  that  the  magnitude  of  the  aerodynamic  force  is  underestimated  by  invoking  the  hyperthermal 
assumption. 


Flat  Plate  -  Lift  Coefficient 


Figure  3.4.  CL  for  Flat  Plate  in  Free-molecular  Flow. 
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Flat  Plate  -  Drag  Coefficient 


Figure  3.5.  CD  for  Flat  Plate  in  Free-molecular  Flow. 


3.2.3  Right  Circular  Cylinder  in  Free-molecular  Flow 

We  now  calculate  the  aerodynamic  force  and  torque  on  the  cylinder  shown  in  Figure  2.7  neglecting 
the  effects  of  the  two  ends  (these  faces  are  simple  flat  plates  which  have  been  treated  previously).  The 
cylinder  has  length  l  and  radius  a  while  its  axis  is  at  an  angle  (3  with  respect  to  the  incident  mean 
molecular  velocity  V.  As  in  the  earlier  hyperthermal  analysis,  it  suffices  to  consider  0</3<  nil .  A 
rectangular  coordinate  system  with  origin  at  the  centroid  O  of  the  cylinder  has  its  z  axis  along  the  axis 
of  the  cylinder  and  x  axis  in  the  plane  determined  by  V  and  the  cylinder  axis.  In  terms  of  the 
cylindrical  coordinates  (r,0,z)  ,  the  lateral  surface  is  described  by  r~a,0<6<2n,-lj2<z^l/2 
and  the  element  of  surface  area  is  given  by  dA  =  ad6dz .  Observing  that 
n  ^  -  cos Oi  sin 0  j  and  v  =  sin/3i  +  cos/3k,  we  have  cosa  =  -sin/Jcos0. 

It  follows  from  eq.(3. 11)  that 

it 

7\=2aljri(ncos6)de  (/rsS sin/3) 

o 

Making  the  substitution  x  =  /r  cos 6 ,  and  introducing  eq.(3.6)  we  obtain 

dx 

Performing  the  integrations,  we  obtain 

%-Fnal  ev/2  [(l  +  )/0  (/u2  /2)  +  f  /,  (/r2  /2)] 

where  /„(  )  denotes  the  modified  Bessel  function  of  the  first  kind  of  order  v. 


2 al  jfe'1’  +*Jnxerf(x) 
4n{  Ju2-x 2 


(3.20) 


In  a  similar  fashion  (introducing  eq.(3.8)),it  can  be  shown  that 


r  /  \  2al  lrX2e  x "+Vrc/2(x  +  2x3)erfx 

%  =  2al  \rJiicos6)cos9ddi  =  -j==—  | - '  v  - dx  i 


Performing  the  integrations,  we  find 


3T  =  ^/p*V/2[(2p2  +3)/0(/x2/2)  +  (2^2  +l)/l(iu2/2)]  i 


Likewise, 


»  t  p  2  /"  1  1 

=  2 a/p/S  f  ri(iucos0)cos26<r/0i  -  — —  f  -T^=^  — *’ +  x  erf  x  Jxi 
J  pSJ0Jp2-x2  W* 


Performing  the  integrations,  we  obtain 


(3.21) 


^  =  ^«//r/5^:/2[(3  +  V)/0(/i72)  +  (V-l)/,(^2/2)] 


Finally, 


=  2aZf r,(pcos0)cos0rZ0i  =  f  ,  ■—  ^dxi  =  —aljX  i 
o  ^  o  -v  ^  ~  x  2 


(3.22) 


(3.23) 


Inserting  the  expressions  from  eqs.(3.20)-(3.23)  into  eq.(3.10)  yields  the  aerodynamic  force  on  the 
lateral  surface  with  components/,  and /  given  by 

/  =  pV2flz||cr„^sinj6-2^(2<7„-a,-4)^sin2)3[(2Ai2+3)/0(^72)  +  (2^+l)/1(/t2/2)]j 
/  =V^pV2a/(7,p'lsin/cos/ev/2[(l  +  /2)/0(yu2/2)  +  /r2/,(/i2/2)]  (3.24) 

Recalling  that  p  =  Ssin/L  we  should  recover  the  hyperthermal  result  derived  earlier  by  taking  the 
limit  of  the  above  expressions  as  p  — >  Utilizing  the  asymptotic  expansions  for  the  Modified 
Bessel  functions 


we  find 


(in2  +3)/„(^V2)  +  (2jJ:  +  l)l,(n,/2)  4 


ljm  (^2  +  1)7o (/*72)  +  /*2/i(a*72)  =  _2_ 

iip*  Vtt 
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Thus  the  limiting  value  of  the  forces  in  eq.(3.24)  as  S  00  do  indeed  degenerate  to  those  given  in 
eq.  (2.16). 

The  axial  and  transverse  forces  (fz/x)  on  the  lateral  surface  of  the  cylinder  are  plotted  in  Fig.  3.6  as  a 
function  of  the  angle  of  attack  (for  a„=a,= 0.7  and  Vw/V=0.05)  with  the  solid  curves  corresponding  to 
eq.(3.24)  with  S=3  and  the  dashed  curves  indicating  the  corresponding  hyperthermal  result  as  given 
by  eq.(2.16).  In  both  cases,  fx  vanishes  at  fi  =  0  and  attains  its  maximum  value  at  P  =  n/2.  We  note 
however  that  at  /3  =  0,  f.  is  non  zero  for  a  finite  speed  ratio  (cr,  *0)  while  it  is  zero  in  the 
hyperthermal  limit.  A  similar  result  was  observed  in  the  previous  case  of  the  flat  plate.  Again,  the 
plots  show  that  the  magnitudes  of  the  aerodynamic  forces  are  underestimated  by  invoking  the 
hyperthermal  assumption. 


Forces  on  Cylinder 


Figure  3.6.  Forces  on  Lateral  Surface  of  Cylinder  in  Free-molecular  Flow 

To  calculate  the  aerodynamic  torque,  we  need  to  evaluate  the  four  integrals  in  eq.(3.13).  It  is  readily 
shown  that  the  integrals  QlxQ^,  and  Q4  all  vanish.  The  remaining  integral  is  simply  related  to  the  above 
value  of  T4. 

g  =  -aTi=-^a2lM\  (3.25) 

The  aerodynamic  moment  about  the  centroid  of  the  cylinder  then  follows  from  eq.(3.12) 

g  =  —  a2lpV2o,  sin  2/S  j  (3.26) 

4 

Interestingly,  the  moment  unlike  the  force,  is  independent  of  the  molecular  speed  ratio  and  of  course 
agrees  with  the  hyperthermal  result  given  in  eq.(2.I7). 
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4.  Rotating  Axisymmetric  Bodies  in  Hyperthermal  Free-molecular  Flow 

So  far  in  our  investigation,  the  bodies  under  consideration  were  not  rotating  with  respect  to  the 
oncoming  (mean)  flow  and  hence  the  surface  relative  velocity  of  the  incident  molecules  did  not  vary 
over  the  impinged  surface  area.  We  now  consider  the  general  class  of  axisymmetric  convex  bodies 
which  are  rotating  about  their  axis  of  symmetry  but  have  arbitrary  orientation  relative  to  the  incident 
molecular  flux.  In  order  not  to  obfuscate  the  effects  of  rotation,  we  restrict  our  attention  to  the  case  of 
a  hyperthermal  free-molecular  flow;  this  assumption  is  removed  in  the  subsequent  section.  General 
expressions  in  terms  of  definite  integrals  are  obtained  for  the  aerodynamic  force  and  moment  resolved 
along  a  set  of  space  fixed  axes.  The  quadratures  are  carried  through  explicitly  for  some  simple 
geometries  and  force  and  moment  components  which  arise  from  the  body's  rotation  are  plotted 
against  the  incident  angle  of  attack. 


4.1  Surface  Geometry 

Figure  4.1  shows  a  body  rotating  about  its  axis  of  symmetry  at  rate  to  relative  to  the  atmosphere.  We 
define  a  spatially-fixed  coordinate  system  {x.y,z}  with  corresponding  unit  vectors  i,  j,  and  k  taking 
the  z  axis  along  the  axis  of  symmetry.  The  incident  molecules  have  velocity  V0  relative  to  the  origin 
O  where  V0  is  at  an  angle  (3  with  respect  to  the  -z  axis  and  is  taken  (without  loss  of  generality)  to  lie  in 
the  y-z  plane.  At  the  element  of  surface  area  clA  instantaneously  at  point  r,  we  construct  a  body  fixed 
frame  with  unit  vectors {t,n,b}  along  the  local  tangent,  normal  and  binormal  directions  i.e.  t  is  along 
fflxr ,  n  points  in  the  direction  of  the  inner  normal  to  the  surface,  and  b  =  t  x  n.  It  proves  convenient 
to  introduce  cylindrical  coordinates  { p,(p,z }  related  to  the  rectangular  coordinates  by 

x  =  ps'uup,  y  =  -pcos<p  (4.1) 

Note  that  since  the  flow  is  symmetric  with  respect  to  the  y-z  plane,  this  choice  of  (p  will  render  many 
expressions  to  be  even  functions  of  (p  thus  simplifying  the  integration.  Since  the  z  axis  is  an  axis  of 
symmetry,  p  will  be  a  prescribed  function  of  z  on  the  surface.  The  position  vector  r  to  a  point  on  the 
surface  is  given  by  the  expression 

r(z,^)  =  psin<pi-pcos<pj  +  zk  (4.2) 

The  unit  vectors  t  and  b  are  along  the  directions  of  dr/d(p  and  dr/dz  respectively,  while  n  =  b  x  t.  If 
we  define  0  as  the  angle  between  the  outer  normal  to  the  surface  and  the  z  axis,  we  can  write 

t  =  cos<pi  +  sinq?  j 

b  =  -  cos  0  sin  (oi  + cost?  cos  <p  j  +  sin0k  (4.3) 

n  =  -sin<psin0i  +  cos(osin0  j- cos  0  k 
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Figure  4.1.  Axisymmetric  Body  in  Free-molecular  Flow 


Thus  the  transformation  matrix  T  from  the  {n.b,t}  to  the  {i.j.k }  basis  can  be  expressed  as 


T  = 


'-sin©sin<p 
sin  ©cos <p 
,  -cos  © 


-cos©sin  (p 

COS©COS<j9 

sin© 


cos  (p' 
sin<p 


(4.4) 


The  element  of  surface  area  is  given  by  dA  =\\dr/d(pxdr/dzd<pdz\\  which  can  be  expressed  in  the 
form 


dA  =  pcscd  d(pdz 


(4.5) 


The  velocity  V  of  the  incident  molecules  relative  to  the  surface  is  given  by 

V  =  V0  -fl)xr 


In  contrast  to  the  non-rotating  case,  the  surface  relative  velocity  now  varies  over  the  surface. 
Expressed  in  the  surface-fixed  triad,  we  obtain 


V  =  V0  (cos  (3  cos  6  +  sin  /3  cos  (p  sin  0)  n  +  V0  (sin  (3  cos  9  cos  (p  -  cos  jS  sin  6)  b 
+  (  Vq  sin/)sin<p-pm)t 


(4.6) 


To  a  high  degree  of  approximation,  we  can  replace  co  with  the  body's  inertial  angular  velocity.  As 
discussed  in  Sec.  2.2,  it  is  important  to  limit  the  aerodynamic  pressure  to  those  portions  of  the  surface 
which  actually  are  impinged  by  the  atmospheric  molecuules.  This  requirement  may  be  expressed  by 
the  inequality  V  n>0.  To  effect  this  limitation,  we  will  insert  the  factor  H(r\)  into  the  surface 
integrals  of  the  aerodynamic  forces  and  torques  to  be  derived  below  where 

?7  =  cos  /3  cos  6  +  cos  (p  sin  (3  sin  6  (4.7) 


4.2  Aerodynamic  Forces  &  Moments 

As  shown  in  Sec.  2.2,  the  incident  flow  has  an  associated  mass  flux  of  pV-n  and  momentum  flux  of 
p(V-n)V.  Denoting  by  pni, pbi,p,i  the  components  of  the  incident  momentum  flux  along  the  local 
normal,  binormal  and  tangent  directions,  we  obtain 

P„,=PV^f 

Pm  =  pVo2?7(cos0cos<psin/3-cos/)sin0)  (4.8) 

P„  =  P  v0  T1  ( V0  sin  (p  sin  j3  -  pco) 

In  addition  to  the  above  incident  momentum  flux,  we  require  an  expression  for  pw  -  the  normal 
component  of  momentum  flux  which  would  be  carried  away  from  the  surface  by  diffusely  reflected 
molecules  in  thermal  equilibrium  with  the  surface.  Utilizing  the  above  expression  for  the  mass  flux, 
we  obtain 


P,=PV0WW  (4-9) 

where  l/w  is  given  by  eq.(2.3).  From  the  definitions  of  the  momentum  accommodation  coefficients 
(eq.(2.1)),  we  can  calculate  the  reflected  momentum  flux  in  terms  of  the  incident  momentum  flux, 
pw,  an  andcr,.  Subtracting  the  reflected  from  the  incident  momentum  flux  yields  the  force  per  unit 
area,  P  ,  exerted  by  the  molecules  upon  the  surface. 

P  =  P  n  +  Phb  +  P  t 

/  (4.10) 

where  PK  =  (2  -  a„  )pni  +  (J„pn„  P„  =  alPhi ,  P,  =  a,p. 
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With  the  aid  of  eq.(4.3).  we  can  now  calculate  the  net  force  F  on  the  body  resolved  along  the  axes  of 
the  fixed  {x.y,z}  coordinate  system . 

=$[^cos^-^,sin£?sin<jc>-  Ph  cos  9  sin  cp]H(ii)d A 
Fy  =  $[^s'n<P  +  P„  s\ndcos(p  +  Phcos9cos(p]H(ri)dA  (4.11) 

F  =$[Phsm9  ~  Pn  cos 9] H(ri)dA 

In  a  similar  fashion,  the  aerodynamic  moment  M  about  the  { }  origin  can  be  calculated  with  the 
aid  of  eqs.(4.2),  (4.3)&  (4.10) 

Mv  =  $[cos  (p  (p  cos  9  -  z  sin  9)  Pn  -  cos  (p(z  cos  9  +  p  sin  9)  Ph  -  z  sin  (p  P,  ]  H{r\)  dA 
My"§[{Pcosd-^ne)^(pP„-{zcos9  +  ps\n9)s\n(pPh+zcos(pP,]H(ri)dA  (4.12) 

M:=§pP,H(n)dA 

The  force  and  moment  expressions  given  above  are  in  a  form  suitable  for  numerical  implementation. 

We  now  proceed  to  expand  the  above  force  and  moment  expressions  so  as  to  reveal  in  an  explicit 
manner  the  effects  of  body  rotation  upon  these  aerodynamic  disturbances.  Due  to  symmetry,  many 
simplifications  will  occur  in  the  process  of  generating  the  final  formulas  thus  allowing  for  closed 
form  solutions  to  several  geometric  primitives. 

Combining  eqs.(4.8)-(4.10)  we  obtain 

p»  =(2-°„)PK2r  +  ex  „pV0Vu,ri 

Ph  =  cr(pV^2?7  (cos0cosq9sin/?-cos/?sin£>)  (4.13) 

P,  =  (T,pV'077(V0  sin  (p  sin  (3-pco) 

We  make  the  following  observations 

1) p  and0  are  independent  of  tp 

2)  q  is  an  even  function  of  <p 

3)  P„  and  Pb  are  even  functions  of  cp 

4)  P,  is  the  sum  of  an  even  and  odd  function  of  tp 

The  region  of  integration  is  described  by:  {-n  <  (p  <  n\  z,  <  z  <  z2}and  the  element  of  surface  area  is 
given  by  eq.(4.5).  In  light  of  the  above  observations,  we  can  express  the  force  components  in  the  form 
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Fx  =  -2a,pVn(o\  J  p2  esc  0  cos  <p  t]  H(ij)  dcp  dz 

7,  0 

h  n 

Fy  =  2PVo  J  |  [sin  0  cos  (p  ((2  -  cj„  )  ?7  +  ct„  VJV0 )  + 

z,  0 

c,  cos  0  cos  cp  (sin  ft  cos  6  cos (p  -  cos  /3  sin  0)  +  cr,  sin  ft  sin2<p  ] 
tj  H(rj)p  esc  6  dcp  dz 


(4.14) 


^2  « 

F  =  2 pV^2  J  J  [crf  sin  0  (sin  /?  cos  <p  cos0  -  cos  ft  sin  0)  -  cos 0  ((2  -  cr ,  )?7  +  an  Vn, /  V0 )] 
J,0_ 

11  H{r\)p  esc  Qdcpdz 


In  a  similar  fashion,  the  components  of  the  moment  can  be  reduced  to 

z-i  n 

=2Pyoj J[cos<p(pcos0  -zsin0)((2- an)r]  +  or,,  VJV0)- 

c,  0 

<7,  cos  (p  (z  cos  0  +  p  sin  0) (cos 0 cos (p  sin  -  cos /? sin  0)  - 
cr,z  sin  /3  sin2<p  ]?j  H{r])pcscG  dcp  dz 

(4.15) 

z2 

Mv  =  -2(7,pV0ffl|Jzp2  cos  <p  ii  H(ii)csc  6d(pdz 

z,  0 
c2  n 

M.  =  -2a,pV0co  j  j  p^r/  H(ri)csc0  dcpdz 

zt  0 

We  observe  that  the  body’s  rotation  has  no  affect  on  Fv.  F,  and  M x  but  introduces  additional 
disturbances  F„  M,  and  M,  proportional  to  the  body’s  angular  rate.  It  is  interesting  to  note  that  the  lift 
force,  F„  acts  in  opposite  direction  to  that  which  would  obtain  under  a  continuum  flow  (Magnus 
effect)  i.e.  we  can  show  that  the  lift  force  corresponding  to  a  free-molecular  flow  acts  in  the  direction 
of  (oxV0.  Taking  co  >  0.  we  will  thus  show  that  Fx  <0.  Based  upon  eq.(4.14),  we  must  prove  that 

*■2  n 

J  p2  esc  0  J  cos  (p  ii  H(ii)d(p  dz  >  0 

z ,  0 

Since  0,/3e(O,;r),  sin/3>0  and  sin0>O.  It  is  therefore  sufficient  to  demonstrate  that  the  inner 
integral  is  nonnegative.  From  eq.(4.7),  we  can  write  r]  =  acoscp  +  c;  a  =  sin/3sin0 > 0,  and 
c  =  cos  p  cos  0 .  It  is  easy  to  show  that 

n  nil 

J  cos  (p  11 H ( rj )  dcp  =  |  [(a  cos  (p  +  c)  H(a  cos  (p  +  c)  -  (- a  cos  (p  +  c)  H(-a  cos  (p  +  c)]  cos  (p  dcp 
0  0 

For  0<cp  <n/2,  acos<p +  c> -a  cos  <p  +  c;  and  of  course  H{  )  is  a  monotonically  increasing  function 

n 

of  its  argument.  It  follows  that  j  cos  cpii  H(ii)  dcp  >0  ( QED ) 

0 
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4.3  Closed  Form  Solutions 


We  provide  below  closed  form  expressions  for  the  aerodynamic  forces  and  moments  on  a  circular 
disk,  sphere,  right  circular  cylinder  and  right  circular  cone  where  each  body  rotates  about  its  axis  of 
symmetry.  Since  results  for  the  non-rotating  case  have  been  provided  in  Sec.  2.3,  we  shall  only 
consider  disturbances  arising  from  body  rotation. 


4.3.1  Spinning  Disk  in  Hyperthermal  Free-molecular  Flow 

Consider  the  disk  x 2  +  v2  <  R2  lying  in  the  plane  z=0.  Restricting  the  angle  of  attack  to  0  <  (3  <  n/ 2 , 
we  have  impingement  only  on  the  top  surface  (0  =  0)  where  r]-cos/3  (see  eq.(4.7)).  The  element  of 
surface  area  is  given  by  dA  =  pdpd(p.  Due  to  the  degenerate  nature  of  this  surface,  we  need  to 
modify  the  form  of  the  integrals  given  in  eqs.(4.14)  &.  (4.15).  The  region  of  integration  will  be 
{0<  p  <  /?;  0  <  (p  <  2?r}  (recall  that  advantage  was  taken  of  symmetry  in  eqs.  (4.14)&;(4.15)  thus  the 
factor  of  2  and  integration  over  0  <  (p  <  n). 

In  R 

Fx(a>)  =  -cr,  pV0(a  cos  ft  j  j  p2  cosy  dp  d(p  =  0 

0  0 

My((o)  =  0 

In  R 

M.(co)  =  -cr,pV'0<y  cos/3  J  J p2dpd(p  = - a,pV0R4co cos/J 

oo  ^ 

Thus  rotation  has  no  affect  upon  the  force  and  produces  a  torque  opposing  the  direction  of  spin. 


4.3.2  Spinning  Sphere  in  Hyperthermal  Free-molecular  Flow 

Figure  4.2  shows  a  sphere  of  radius  R  rotating  about  the  z  axis  with  angular  rate  (0.  The 
hemispherical  portion  of  the  sphere  impinged  by  the  oncoming  flow  (incident  at  angle  p  with  respect 
to  the  z  axis)  is  not  conveniently  described  by  the  set  of  cylindrical  coordinates  { p,(p,z }  defined  in 
Figure  4.1.  We  therefore  transform  to  a  new  rectangular  coordinate  system {X,  Y,Zj  by  rotating  about 
the  x  axis  by  angle  p.  A  set  of  spherical  coordinates  {/?,©, <£>}at  the  center  of  the  sphere  is  defined  by 
the  transformations 


X  =  RsinQsinO,  Y  =  -/?sin0cos<t>,  Z=R  cos©  (4.16) 

The  element  of  surface  area  is  given  by  dA  =  R2  sin 0c/©r/O  and  the  region  of  integration  (impinged 
portion  of  surface)  is  conveniently  described  by  {-n  <  O  <  n\  0  <0  <71/2].  In  order  to  apply  the 
general  formulas  given  by  eqs.(4.14)&(4.15),  we  need  to  relate  the  cylindrical  coordinates  {p,(p,z}  to 
the  spherical  coordinates  {/?,©,0}.  The  transformation  between  the  two  sets  of  rectangular 
coordinates  is  given  by 


(x, y,z)  =  (X, cos/3 T-sin/JZ,  sin/JT  +  cos/3Z) 
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Figure  4.2.  Rotating  Sphere  in  Free-molecular  Flow 


Introducing  eqs.(4.1)  &  (4.16)  into  the  above,  we  obtain  the  relations 

psintp  =  ft  sin  0  sin  O 

pcos<p  =  ft(sin)3cos0  +  cos/3sin0cos<E>)  (4.17) 

z  -  ft(cos/3cos0~sin/lsin0cosd>) 

Recalling  the  definition  of  0,  we  see  that  for  the  present  case  of  the  sphere,  cos d  =  z/R  and 
sin 6  =  p/R.  Utilizing  these  results  in  eq.(4.7),  we  obtain,  as  expected, 77  =  cos 0.  The  flow  being 
symmetric  with  respect  to  the  y-z  plane  implies  that  all  integrands  in  eqs.(4.14)  &  (4.15)  are  even 
functions  of  O.  Performing  the  integrations,  we  obtain. 

Fx  (co)  =  -— a,pV0R3co  sin  j6 

My(co)  =  -—a,  pV0R4(O  sin  2yS  (4.18) 

8 

Mz(co)  =  -—a,  pV0R4co  (5  -cos  2)3) 

8 

The  two  torque  components  (scaled  by  no ,pV0R4 a)  /  8)  are  plotted  below  as  functions  of  the  angle  of 
attack.  Note  that  the  rotation  of  the  sphere  results  in  a  torque  opposing  the  spin  direction  over  the 
entire  range  of  angle  of  attack. 
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Torques  on  Rotating  Sphere 


Figure  4.3.  Spin  Moments  on  Sphere  in  Hyperthermal  Free-molecular  Flow 


4.3.3  Spinning  Cylinder  in  Hyperthermal  Free-molecular  Flow 

With  reference  to  Fig.  4.1.  a  right  circular  cylinder  of  radius  R  and  length  L  has  its  centroid  at  the 
origin  of  the  Cartesian  coordinate  system  {x,y,z}.  The  cylinder  is  rotating  about  its  axis  ( z )  at  angular 

rate  0)  with  an  impinging  free-molecular  flow  incident  at  angle  of  attack  |3.  On  the  lateral  surface  of 
the  cylinder.  0  =  7t/2,  r]  =  sin)3cos<p  and  for  0 <(3<n,  77 >0  if  |<p|<7r/2.  It  follows  from  eqs.(4. 14) 
&  (4.15)  that  the  spin  induces  a  force  on  this  portion  of  the  surface  given  by 

Flx  =  -^-cocr,  pV0R2L  sin 

and  a  moment  about  the  center  of  the  cylinder 

M'y=  0 

Ml  =-2£O(T,pV0/?3Lsinj8 

If  we  restrict  0<  /?<7t/2,  then  there  will  be  no  impingement  on  the  bottom  surface  (z  =  -L/2)of  the 
cylinder.  The  force  and  moment  on  the  top  portion  of  the  cylinder  follow  from  the  analysis  above  of 
the  spinning  disk.  Thus  the  net  force  and  moment  on  the  cylinder  due  to  its  rotation  are  given  by 
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(4.19) 


Fx(co)  =  -y  oxj^VqR2  Lsm  /3 
My(co)  =  0 

Af.(cu)  =  --i(«cr,pF0/?3  (7z:/?cos/3-!-4Lsin/3) 

The  z  component  of  the  moment  (scaled  by  -coa,pV0R3 /2)  is  plotted  below  as  a  function  of  the 
angle  of  attack  for  i)  R-l,L=5  and  ii)  R=5,L=l. 


-m.  Torque  on  Rotating  Cylinder 


Figure  4.4.  Spin  Moments  on  Cylinder  in  Hyperthermal  Free-molecular  Flow 


4.3.4  Spinning  Cone  in  Hyperthermal  Free-molecular  Flow 

Figure  4.5  shows  a  right  circular  cone  of  base  radius  R  and  height  L  with  the  origin  of  the  Cartesian 
coordinate  system  {x,y,z}  at  the  center  of  its  base.  The  cone  is  rotating  about  its  axis  (z)  at  angular 

rate  co  with  an  incident  free  molecular  flow  V0  at  angle  of  attack  (3.  We  restrict  (3  to  the  range  (0,zr/2) 
so  that  there  will  be  no  impingement  on  the  base  of  the  cone.  The  angle  8  between  the  outer  normal  to 
the  surface  and  the  z  axis  is  given  by  0  =  7r/2-y  where  y  is  the  half-angle  of  the  cone  i.e. 
tan y-R/L.  The  cylindrical  coordinate  p  is  given  by  p=(L-z)tany,  and  from  eq.(4.7), 
7]  =  cos ^3 sin y  +  sin  jScosy  cos<p.  To  determine  when  77  is  positive,  it  proves  convenient  to  consider 
the  two  cases  1)  0  <  (5  <  y  and  2)  y  <  n/2. 
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Figure  4.5.  Rotating  Cone  in  Free-molecular  Flow 


•  Case  1 .  0  <  P  <  y 

Here  we  will  have  sin  /3  <  sin  y  and  cos/3>cosy  so  that  cosjSsiny  >  cosy  sin  /3  >  0  and  hence  ?7  >  0 
for  all  (p.  Performing  the  integration  in  eqs.(4.14)  &  (4.15)  we  obtain 

7t 

Fx(a>)  ~-—a,pV0 R3  cot  y  0)  sin  ft 

My(co)  =  jR cot y  Fs((o)  (4.20) 

M,  (<y )  =  -  ^  a,  pV0  R*(o  cos  p 


•  Case  2.  y  <) 3<  n/2 

In  this  case,  cos (3<  cosy  and  siny  <sin/S  therefore  77  is  positive  at  <p  =  0  and  negative  at  (p  =  7t.  It 
can  be  shown  that  r\[(p)  has  a  single  root  on  [0,7r]  at  (p  =  (p *  where 


< P 


_  f 

=  cos'1 

V 


tanyN 

tan/?. 
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Hence  the  integration  limits  are  now  {0<z<L;0<<p<aw.  Since  lim  (p*  =n,  we  will  have  a 
continuous  transition  across  the  point  (3=y.  Performing  the  integration  in  eqs.  (4.14)&(4.15)  we  obtain 

Fx(co)--^alpV0R3 co(<p*  sin/3cot y  +  cos/3sin (p  ^ 

My(co)  =  -^coty  Fx{co)  (4.21) 

Mz (a>)  =  — ^ atpV0 RA  co[(p*  cos  +  cot  y  sin  /3 sin  (p' ) 


Fig.  4.6  is  a  plot  of  the  spin  induced  force  Fx.(a>)  (scaled  by  -a,pV0R3col  3)  as  a  function  of  the 
angle  of  attack  for  a  cone  with  half-angle  of  30°.  It  is  interesting  to  note  that  the  magnitude  of  this 
force  starts  out  as  an  increasing  function  of  the  angle  of  attack  (P)  but  does  not  assume  its  maximum 
value  at  /3  -  90° . 

The  spin  induced  torques  My(co)  and  (scaled  by  -cy^VqR4 co)  are  shown  in  Fig.  4.7  as 

functions  of  the  angle  of  attack  on  the  30°  cone.  We  note  that  the  magnitude  of  Mz  is  significantly 
larger  than  that  of  My . 
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Figure  4.6.  Spin  Force  on  Cone  in  Hyperthermal  Free-molecular  Flow 


Figure  4.7.  Spin  Moments  on  Cone  in  Hyperthermal  Free-molecular  Flow 


5.  Rotating  Axisvmmetric  Bodies  in  Free-molecular  Flow 

We  proceed  to  generalize  the  analysis  presented  in  the  previous  section  for  a  body  rotating  about  an 
axis  of  symmetry  in  a  free  molecular  flow  by  removal  of  the  hyperthermal  assumption.  General 
expressions  in  the  form  of  surface  integrals  are  obtained  for  the  aerodynamic  forces  and  moments 
which  reduce  to  the  previously  obtained  hyperthermal  results  in  the  case  of  an  infinite  molecular 
speed  ratio.  Again  it  is  found  that  the  body  rotation  gives  rise  to  additional  force  and  moment 
components  proportional  to  the  body’s  angular  rate.  Although  quite  intricate,  closed  form  solutions 
are  obtained  for  some  simple  body  shapes  and  errors  introduced  by  the  hyperthermal  assumption  are 
investigated  by  examining  the  variations  in  the  aerodynamic  disturbances  over  a  range  of  molecular 
speed  ratios  and  angles  of  attack.. 


5.1  Mass  and  Momentum  Flux 

As  depicted  in  Figure  4.1.  an  axisymmetric  body  rotates  at  angular  rate  to  about  its  axis  of  symmetry 
relative  to  the  atmosphere  We  now  suppose  that  a  random  velocity  u  described  by  the  Maxwell- 
Boltzmann  distribution  is  superimposed  on  the  mean  atmospheric  flow.  Define  a  spatially-fixed 
coordinate  system  {x,y,z}  with  corresponding  unit  vectors  i,  j,  k  such  that  the  z  axis  is  along  the  axis 
of  symmetry.  The  incident  molecules  have  a  mean  velocity  V0  relative  to  the  origin  O  where  V0  is  at 
an  angle  (3  with  respect  to  the  -z  axis  and  is  taken  (without  loss  of  generality)  to  lie  in  the  y-z  plane. 
The  velocity  V  defined  by  V  =  V0  -  tax  r,  is  the  velocity  of  the  incident  molecules  relative  to  the 
surface  at  position  r  excluding  the  random  thermal  motion  of  the  atmosphere.  The  components  of  V 
along  the  local  {n,b,t}coordinate  directions  (see  Sec.  4.1)  are  given  respectively  by 


Vn  =  V0  (cos  ft  cos  6  +  sin  /?  sin  9  cos  tp) 

Vb  =  Vo(sin/3cos0cos<p-cos/2sin0)  (5.1) 

V,  =  V0  sin/lsin (p-a>p 


The  number  of  molecules  per  unit  volume  with  random  velocities  between  u  and  u+du  is  given  by 
eq.(3.1) 


mV„n 


3_3/2 


exp 


.2  3 


dut  diin  du 3 


J 


where  p  is  the  gas  (atmospheric)  density,  m  is  the  mass  of  a  single  gas  molecule,  and  V„  is  the  most 
probable  speed  of  the  gas  molecules.  The  random  velocity  u  is  with  respect  to  the  mean  motion  of  the 
gas  and  hence  the  velocity  v  of  the  incident  molecules  relative  to  the  surface  including  the  random 
thermal  motion  is  given  by  v=V+u.  Resolving  velocities  along  the  local  {n,b,t}coordinate  directions, 
it  follows  that  the  number  of  molecules  per  unit  volume  with  surface  relative  velocity  between  v  and 
v+dv  is  given  by  f(vn,vb,vl)dvndvbdvl  where 


f{Vn’VbV,)  =  - 


-exp 


V  ■  V 

v: 


exp 


2V -v-v-v 


(5-2) 


Following  the  argument  given  in  Sec.  3.1,  the  net  mass  flux  8Q  (mass  delivered  per  unit  time  per  unit 
surface  area)  from  all  incident  molecules  is  given  by 
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SQ  =  mj  jjvn  f{Vn  ’  Vfc  ’  Vf  )dv„  dvb  dvt 


which  can  be  expressed  in  the  form 


8Q  =  -^T™  P 


v-v 

v '  K2  j 


2V,v , -v,2^  7 
exp  . '•  Idv,  j 


exp 


(  2  Vbv„-vl' 
V2 

\  Va 


dvb\Vn 


exp 


^Vv-v2)  , 

-  n  ",  "  </v„ 

v  K2  J 


Performing  the  integration,  we  obtain 


6Q  =  pVori(57?)/S  (5.3) 

where  S  =  V0/Va  (molecular  speed  ratio;  S  =  °°  in  the  hyperthermal  limit),  T,  (jc)  is  given  by  eq.(3.6), 
and  t|  is  defined  in  eq.(4.7).  The  incident  momentum  flux  p,  due  to  all  impinging  molecules  is  given 
by  the  expression 


v„n  +  vfib  +  v,t)v„  /(v„  ,vb,v,  )dvn  dvb  dv, 


Utilizing  eq.(5.2)  and  denoting  the  components  of  p;  along  the{n,b,t}  coordinate  directions  by 
^respectively,  we  find 


p*=« 4?exp 


( 


V-V^l?  (2Vbvb-vl' 


\  r« 


J 


exp 


\dvk 


\  va 


J 


Jv>  p| 

f  2Vv  -v„2^l 

V2 

0 

\  Va  J 

dv„ 


Pbi  = 


P 


P-,=^Tex  P 


^eXP(“^"]  I  VfceXP^~y2  ^  j^Vi>  J  eXp[2V,y2  V'  jJvJV„eXP 

V-V'lf  {2Vbvb-v2b)J  }  (2V,v,-Vj)J  7 

exp[ — W~^)dVb  J  v,exp[  v2  jdv,J  v"exp 


2V  v  -  v 

n  n _ n 


2  'N 


dv. 


\  V«  J 


2V„v  -v  V 

"  w2  dv„ 
Va 


Performing  the  integrations  and  simplifying,  we  obtain  the  following  expressions  for  the  components 
of  the  incident  momentum  flux. 


Pm  =  PVo  r2(Sr])/ S2 
Pbi  =  pV0Vbri{Sr1)/S 
Pu  =  PVq  V,  r^Srfj/S 


(5-4) 


where  r2(jc)  is  given  by  eq.(3.8). 


We  also  require  an  expression  for  pw  -  the  normal  component  of  momentum  flux  which  would  be 
carried  away  from  the  surface  by  diffusely  reflected  molecules  if  they  were  in  thermal  equilibrium 
with  the  surface.  We  can  write  pw  =  VwSQ  where  Vw,  given  by  eq.(2.2),  is  average  normal  velocity 
for  this  class  of  molecules.  Invoking  eq.(5.3),  we  obtain 

Pw  =  PVqVw  r^Srf/S  (5.5) 
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5.2  Aerodynamic  Forces  &  Moments 

The  force  per  unit  area,  P,  exerted  by  the  molecules  upon  the  surface  is  given  by  eq.(4.10)  in 
conjunction  with  eqs.(5.4)  &  (5.5)  above. 


p„=(  2  -  H2  r 2(sv)/s2  +  anpv0vw  r,  (srf/s 

Pb=(J,pV0Vbri(Sri)/S  (5.6) 

P<=°,pv0v,r](sii)/s 

The  expressions  for  the  net  aerodynamic  force  and  torque  along  the  spatial  axes  in  the  present 
situation  is  given  by  eqs.(4.11)  &  (4.12)  with  H(ji)  replaced  by  unity  and  are  in  a  form  suitable  for 
numerical  evaluation.  We  now  proceed  to  expand  these  force  and  moment  expressions  so  as  to  reveal 
in  an  explicit  manner  the  effects  of  body  rotation  upon  these  aerodynamic  disturbances.  Considerable 
simplification  in  evaluation  of  the  integrals  can  be  realized  by  making  the  following  observations. 

1)  p  and  0  are  independent  of  (p 

2)  Pn,ar\dPb  are  even  functions  of  (p 

3)  Pt  is  the  sum  of  an  odd  and  even  function  of  (p 

The  region  of  integration  is  described  by:  \-n  <  <j o<n\zx  <z<  z2}  and  the  element  of  surface  area  is 
given  by  eq.(4.5).  In  light  of  the  above  observations,  we  can  express  the  force  components  in  the  form 


-2  7* 

Fx  =  ~2 C7,pV0  co/S  J  Jp2  cscO cos cpT^Si^dcpdz 

t,  o 

*2  n 

Fy  =  2pV^2  J"  JjJtf,  sin /3 sin2  ip +  cr,^ cos 0 cos ^)  +  cr„  —  sin0cos<pjr,(S?7)/S- 
=.  o  0 

(2  -  an )  sin  6  cos  (pF2(Sri)/ S2 1  p  esc  6  depdz 


*2  n 

F.  =  2pVg  J*  J[lCT  ,B,  sin  9  -  an  — -  cos  o)  T,  ( Sri)/S  -  (2  -  <Jn )  cos  6  T2  (Sp)/S2  lpcsc  9  depdz 
„  kn  ' 


where  £  =  sin  /3  cos  0  cos  (p  -  cos  /3  sin  6 


The  components  of  the  aerodynamic  moment  about  the  origin  of  the  spatial  axes  are  given  by 


Mx  =  2pVj  Jf{h  Vw /  V0  (p  cos  6  -  z  sin  9)  cos  (p  -  a,  [z  cos  9  +  p  sin  £>)£ cos  (p 

c,  o 

-0‘(jsinj3sin2^]r,(5??)/S  + 

(2-°»)(pcos0-  z  sin  0)  cos  (p  r2  (S  i])/ S 2  }p  esc  6  d(pdz 

?r  (5-8) 

2  eseQ  co&(pT^Sri)d(pdz 

*1  o 
22  « 

M,  =  -IO'PVqCo/S  j*  J  p3  esc  9T\(Sii)d(pdz 

M  o 

As  in  the  hyperthermal  case,  the  body’s  rotation  has  no  affect  on  F, F,  and  MA  but  introduces 
additional  disturbances  Fx,Mv  and  M.  proportional  to  the  rotational  rate.  Additionally,  we  can  again 
show  that  the  lift  force  acts  in  the  direction  of  <axV0;  which  according  to  Fig.  4.1,  requires  Fx  <0 
forty  >  0.  Hence  we  must  prove  that  (see  eq.(5.7)) 

Z2  n 

J  p2  esc  6  J  cos  (p  T,  ( Sr\)dcp  dz>  0 

c,  o 

Since  0,/3e(O ,n).  sin0>O  and  sin/?>0.  It  is  therefore  sufficient  to  demonstrate  that  the  inner 
integral  is  nonnegative.  From  eq.(4.7)  we  can  write  r}  =  acos,(p  +  b\  where  a  =  sin/3sin0> 0,  and 
b  =  cos  / 3  cos 6 .  It  is  easy  to  show  that 

n  nil 

Jcospr^S??)^  J  [r,(SZ>  +  Sa  cos  <p)  -  F,  (Sb  -  So  cos  <p)]cos^>  d(p 
0  0 

Now  the  function  r,(A')  is  an  increasing  function  of  x  and  Sb  +  Sa  cos  (p  >  Sb  -  Sa  cos  (p  for 

it 

0<(p<n/2.  It  follows  that  Jcos^r,(5?7)t/^  >  0  ( QED ) 

o 


Recalling  the  two  limits  given  by  eqs.(3.14)  &  (3.15) 


lim  Lipl  =  n  //( 77),  Jim 


=  rf  H(v) 


where  H(x) is  the  Heaviside  step  function,  it  is  easy  to  show  that  the  expressions  (5.7)  &  (5.8)  reduce 
respectively  to  the  hyperthermal  results  obtained  previously  in  eqs.(4.14)  &  (4.15). 
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5.3  Closed  Form  Solutions 

Although  the  calculations  are  rather  intricate,  we  will  find  closed  form  expressions  for  the 
aerodynamic  forces  and  moments  on  a  circular  disk,  sphere,  right  circular  cylinder  and  right  circular 
cone  where  each  body  rotates  about  its  axis  of  symmetry.  Since  results  for  the  non-rotating  case  have 
been  provided  in  Sec.  3.2,  we  shall  only  consider  disturbances  arising  from  body  rotation.  It  will  be 
demonstrated  that  the  resulting  expressions  reduce  to  those  obtained  in  Sec.4.3  (hyperthermal 
assumption)  as  the  molecular  speed  ratio  tends  to  infinity. 

5.3.1  Spinning  Disk  in  Free-molecular  Flow 

We  consider  the  disk  x 2  +  y2  <  R2  lying  in  the  plane  z=0.  Note  that  unlike  the  hyperthermal  case, 
even  if  we  restrict  the  angle  of  attack  to  0  <  ji  <  n/2  the  disk  experiences  impingement  on  both  the 
top  and  bottom  surfaces  due  to  the  random  thermal  motion  of  the  molecules.  On  the  top  surface 
8  =  0,t]  =  cos P  while  on  the  bottom  surface  6  =  n,rj  =  -cos  ft .  Due  to  the  degenerate  nature  of  this 
surface,  we  need  to  modify  the  form  of  the  integrals  given  in  eqs.(5.7)  &  (5.8)  by  replacing  the 
element  of  surface  area  (J5  esc  ddcpdz)by  pdpdcp  and  integrating  over  the  region 
{0  <  p  <  R-,0  <(p<  n}.  We  find  that  both  Fk(cd)  and  MY(a>)  vanish  while  the  disk  rotation  produces  a 
torque  opposing  the  direction  of  spin  given  by 

M.  (a)  =  - 1  a,  p Vq  Ra(o  [r,  (S  cos  p)  +  r,  (-5  cos  0)}/S  (5.9) 

Now  lim  M.  =-n/2a,pV0R4(X)  cos/3  in  agreement  with  the  hyperthermal  result  obtained  in  Sec. 

5— 

4.3.1.  At  the  other  extreme,  when  S- 0  (Va  finite,  V0  =  0)  we  have  M.  -~^/2<7lpVaRi(0. 

Figure  5.1  shows  the  variation  of  M.  (scaled  by  -7K7,pV0R4CL)/2)  with  angle  of  attack  for  the  two 
cases  S  =  2and  S  =  It  is  observed  that  the  hyperthermal  approximation  breaks  down  at  high  angles 
of  attack;  at  ft  =  90°  it  yields  zero  torque  while  a  nonzero  value  obtains  under  a  finite  molecular 
speed  ratio. 

In  Figure  5.2,  we  have  plotted  the  percentage  error  in  M.  corresponding  to  the  hyperthermal 
approximation  as  a  function  of  S  for  several  values  of  angle  of  attack.  Again,  the  error  increases  with 
angle  of  attack  and  is  100%  at  ft  =  90°. 
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Moment  on  Rotating  Disk 


Torque  Error  for  Rotating  Disk 


Figure  5.2.  Disk  Hyperthermal  Spin  Moment  Error 
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5.3.2  Spinning  Sphere  in  Free-molecular  Flow 

Consider  a  rotating  sphere  of  radius  R  in  a  free  molecular  flow  as  shown  in  Fig.4.2.  Following  the 
procedure  in  Sec.  4.3.2,  we  rotate  to  a  Cartesian  frame  {X,  Y,  Z }  so  that  the  mean  flow  velocity  V()  is 
parallel  to  the  (-Z)  axis  and  introduce  a  system  of  spherical  coordinates  {/?,©,$}  defined  by 
eqs.(4.16).  The  flow  is  symmetric  with  respect  to  the  {F,  Z}  plane,  hence  the  integrands  in  eqs.(5.7)  & 
(5.8)  are  even  functions  of  d>.  In  the  integrands,  we  set  17  =  cos©,  p  cscd  dcpdz  — >  R2  sin©r/0r/O 
and  make  use  of  eqs.(4.17)  while  first  performing  the  integration  over  O  (recall  that  due  to  the 
random  thermal  motion,  we  must  integrate  over  the  entire  surface  of  the  sphere). 

The  lift  force  Fy(co)  reduces  to  the  definite  integral 

Fx  =  -2narpV0R2a)S^  sin  ft  j  sin©cos©r,(Scos©)£/© 

Making  the  substitution  x  =  S cos©  and  invoking  eq.(3.16)  we  obtain 

Fx  =  -2 tt/3 cr, pVQ R^co  sin  p  (5.10) 

Interestingly,  the  lift  force  is  independent  of  S  and  of  course  agrees  with  the  hyperthermal  result 
obtained  in  Sec.  4.3.2. 


The  expression  for  the  moment  Mv(tu)  can  be  reduced  to 

Mv  =  -7rcr,pV'0/?4(yS~1  s'mip  J  sin  ©(cos2  0-1/2  sin2  ©jr^S  cos  Q)dQ 
Performing  the  same  substitution  as  above,  we  obtain  the  final  expression 


where 


Mv=-n/Sa,pV0R4(oJl(S)sm2P 


2S2  -3  .  4S4  -4S2  +3 

(5)_2V^53eXP(  5  )+  454  6  (  ) 


(5.11) 


Since  limJ?l(S)  =  0  and  Hm  3l(S)  =  1,  this  component  of  the  aerodynamic  moment  vanishes  when 

the  mean  flow  velocity  is  zero  (with  random  thermal  motion)  and  agrees  with  the  result  obtained 
previously  under  hyperthermal  conditions. 

Finally,  the  moment  Mz{a>)  can  be  reduced  to  the  form 

M,  =-n/4<JlpV0Ri<x>S~i j  [5-cos2j0-(l  +  3cos2j8)cos2©]sin0ri(S'cos0)z?© 

Performing  the  integration,  we  find 


where 
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Since  ||m{j?l(S),S(5)}  =  {1,5},  we  recover  the  hyperthermal  limit  as  obtained  in  Sec.  4.3.2.  To 
obtain  the  limiting  value  of  Mz  as  S  (V^)  approaches  zero,  we  replace  V0  with  SVn  and  note  that 
li^i{-5zl(S),SS(S)}  =  jo,32y^3V^)J.  Hence,  under  a  zero  mean  flow  the  sphere  experiences  a  nonzero 
moment  when  spinning  in  a  Maxwellian  atmosphere. 

YimM.  =-4*Jn/3cr,pV„RiCQ 
o  *  n 


The  functions  ^l(S)  and  S( S)  are  plotted  below  in  Fig.  5.3.  Since  0 <  <  1  for  0<S<°°,  the 

magnitude  of  My  is  overestimated  by  the  hyperthermal  approximation.  Figure  5.4  shows  the  variation 

of  \M\  (scaled  by  n/&(7lpV0R*a))  with  angle  of  attack  for  the  two  cases  S=5  and  S=°o.  In  contrast  to 
the  rotating  disk  (Fig.  5.1),  the  torque  magnitude  increases  with  P  and  the  difference  between  the  two 
cases  are  approximately  uniform  over  the  full  range  of  angle  of  attack..  The  plots  contained  in  Figure 
5.5  show  the  percentage  error  in  M.  based  on  the  hyperthermal  approximation  as  a  function  of  S  for 
selected  values  of  the  angle  of  attack  (P).  As  opposed  to  the  problem  of  the  rotating  disk  (see  Fig. 
5.2),  the  error  is  not  greatly  influenced  by  the  angle  of  attack  and  now  decreases  with  increasing  p. 


A(S)  B (S) 


Figure  5.3.  The  Functions  .Ti(S)  &  S(S) 
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-Mz  Moment  on  Rotating  Sphere 


Figure  5.4.  Spin  Moment  A/,(co)  on  Sphere  in  Free-molecular  Flow 


AM*  Torque  Error  for  Rotating  Sphere 


Figure  5.5.  Sphere  Hyperthermal  Spin  Moment  Error 


(deg. ) 
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5.3.3  Spinning  Cylinder  in  Free-molecular  Flow 

With  reference  to  Fig.  4.1,  a  right  circular  cylinder  of  radius  R  and  length  L  has  its  centroid  at  the 
origin  of  the  Cartesian  coordinate  system  {x,  y,  z} .  The  cylinder  is  rotating  about  its  axis  (z)  at  angular 

rate  to  with  an  impinging  free-molecular  flow  incident  at  angle  of  attack  (3  with  respect  to  the  -z  axis. 
Due  to  the  random  thermal  motion  of  the  molecules,  impingement  occurs  on  both  the  top  and  bottom 
faces  of  the  cylinder  in  addition  to  its  lateral  surface.  We  will  use  the  superscripts  T,  7’  &  ‘ b ’  to 
distinguish  forces/moments  on  the  lateral,  top  and  bottom  surfaces  respectively. 

On  the  lateral  surface  we  have  p  =  R,  8  =  n/2  and  77  =  sin  >3  cos  (p  (see  eq.(4.7)).  Computing  the  lift 
force  from  eq.(5.7)  and  performing  the  z  integration,  we  obtain 

Fl(co)  =  -2a ,pV0R2La>S^  [  cos qo T, (^u cos q?) ^r/qo  (p  =  S sin/3) 

J  0 

The  integral  appearing  in  the  above  expression  has  been  previously  evaluated  (see  eq.(3.23))  and 
leads  to  the  result 

F’((o)  =  -^a)(JlpV0R2L  sin/J 

As  in  the  previous  example  of  the  sphere,  the  lift  force  is  independent  of  S  and  agrees  with  the 
hyperthermal  expression  obtained  in  Sec.  4.3.3. 

Calculating  the  moments  on  the  lateral  surface  from  eqs.(5.8)  we  find  that  vanishes  while 

A/' (to)  can  be  expressed  in  the  form 

M'z(co)  = -2atpV0R2LcoS~'  [  F,(pcos (p)dcp 

J  0 

Appealing  to  eq.(3.20)  for  evaluation  of  the  integral,  we  obtain 

M[{co)  =  -4na,pV0RiLa)S~'  <T/<V:  [(p2 +l)l0(p2 /2 )  +  p2  /,(p2/2)] 

where  /v,(  )  denotes  the  modified  Bessel  function  of  the  first  kind  of  order  V.  Exploiting  the  second 
limit  following  eq.(3.24),  we  recover  the  hyperthermal  result  obtained  in  Sec.4.3.3.  At  the  other 
extreme,  we  easily  can  demonstrate  that 

lim  M'Jco)  =  ~4n  <7,pVriR2Lco 
s->  0  ~  '  ‘ 


On  the  top  and  bottom  surfaces  of  the  cylinder  we  have  6  =  0, 77  =  cos 0  and  6  =  n,rj  =  -cos/3 
respectively.  Adopting  the  results  obtained  in  Sec. 5. 3. 1  for  circular  disk,  we  find: 

On  the  top  surface  F'.  =  M'y  =  0 

Mi  (to)  =  -  n/2  a,pV0R*co  S_lr,  (Scos  0) 

On  the  bottom  surface  =  Mby  =  0 

Mbz((o)  =  -  n/2  (T,pVQR4o)  S_l r,  (-Scos  0) 

Adding  the  contributions  from  all  three  surfaces,  we  obtain  the  net  aerodynamic  force  and  torque 
(about  the  centroid  of  the  cylinder)  due  to  rotation. 
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(5.13) 


Fx(a>)  -  -n/2atpV0co  R2L sin  (3 
My(<o)  =  0 

Mz (to)  =  -c T,pV0R3co  S~]  |y /?[r,  (5cos /3)  +  T,  (-5cos  /?)]  + 
4nLe-»1/2[(p2  +  l)l0(p2/2)  +  p2l , (^2/2)]  } 


For  the  purposes  of  numerical  evaluation  we  take  the  cylinder  dimensions  to  be  R- 1  and  L= 5. 

Figure  5.6  shows  the  variation  of  |M,|  (scaled  by  a,pV0CL))  with  angle  of  attack  for  the  two  cases  5=3 
and  S  =  It  is  seen  that  the  error  in  the  hyperthermal  assumption  is  most  pronounced  at  small  angles 
of  attack.  The  Plots  Contained  in  Figure  5.7  show  the  percentage  error  in  Mz  based  on  the 
hyperthermal  approximation  as  a  function  of  S  for  a  range  of  angles  of  attack  between  0  and  90 
degrees.  Again,  it  is  observed  that  the  errors  are  largest  at  small  angles  of  attack  and  are  non- 
negligible  even  at  5=8  for  small  p. 


-mz  Moment  on  Rotating  Cylinder 
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5.3.4  Spinning  Cone  in  Free-moiecular  Flow 

We  consider  a  right  circular  cone  as  shown  in  Fig.  4.5  in  which  a  free  molecular  flow  with  mean 
velocity  V0  is  incident  at  angle  3  to  its  axis.  The  cone  has  base  radius  R  and  height  L  and  is  rotating  at 
the  angular  rate  (0  about  its  axis.  Due  to  the  random  thermal  motion  of  the  molecules,  impingement 
occurs  on  the  bottom  surface  of  the  cone  in  addition  to  its  lateral  surface  for  any  value  of  p. 

On  the  lateral  surface:  6  =  n/2 -y,  p  =  R(l - z/L),  and  t]  =  cos /3sin y  +  sin  P cos y  cos (p 
where  y  is  the  half-angle  of  the  cone  (tan  y  -  R/L ) 

Computing  the  spin  induced  lift  force  and  moments  from  eqs.(5.7)-(5.8)  and  performing  the  z 
integration,  we  obtain 

F[ (co)  =  - — a t pV0/?3 co esc y  —  f  coscpT^Srjjdcp 
3  o 

M\  (co)  =  j  cot  yFl  (co) 

Mz(co)  =  -^ °,Pvo RAu esc y  \  f * T, (Sr]) dtp 

On  the  bottom  surface  of  the  cone,  the  force  and  moment  calculation  follows  that  performed  in  Sec. 
5.3.1  with  0  =  7rand  77  =  -cos/3.  Adding  this  contribution  to  the  above  yields  the  net  aerodynamic 
force  and  moment  (about  the  center  of  the  base)  due  to  rotation. 

Fx(a>)  =  ~a,pV0l^(OcscY  ~  f  cos<pr,(S7 ])dtp 

My(a> )  =  -^coty  Fx(co)  (5.14) 

Mz(co)  =  - ^f0>  r, (-Scos  P)  +  CSC y  [  r,  (St]) dtp } 

Utilizing  the  limit  in  eq.(3.14),  we  can  show  that  these  expressions  simplify  to  the  hyperthermal 
results  obtained  previously  in  Sec.4.3.4.  In  general,  the  integrals  appearing  in  the  above  formulas 
must  be  evaluated  numerically.  We  can  however  obtain  closed  form  expressions  for  the  special  cases 
of  longitudinal  (P=0)  and  transverse  (P  =  jt/2 )  flow. 

When  P=0,  we  readily  find  that  Fx(co)  =  =  0  and 

Mz(co)  =  -^a,pV0 Ra(0 S~]  [P, (-S)  +  esc y T, ( S sin y )] 

When  P  =  n/2,  we  can  show  that 

Fx(co)  =  -n^G^VQR^acoiy 
My((o)  =  -  n/2A  <J,pVQRA  <x>  cot2y 

Mz(co)  =  -4n  /4a,pV0R4co  5_1|l  +  cscy  exp^-K2/2^\  +  K2yo(K2/2')  +  K2I]  (/c2/2)j| 

where  K  =  Scosy  and/v(  )  denotes  the  modified  Bessel  function  of  the  first  kind  of  order  v.  We 
note  that  for  P=0  and  n/2  both  Fx  and  Mv  are  independent  of  S  while  Mz  is  a  decreasing  function  of 
S. 
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The  plots  in  Figure  5.8  show  the  relative  error  in  Fv(ft))  for  the  hyperthennal  approximation  for  (3=5, 
50,  &  70  degrees.  It  is  observed  that  the  error  decreases  with  increasing  angle  of  attack  for  a  given 
molecular  speed  ratio.  This  is  consistent  with  the  fact  established  above  that  Fx  is  independent  of  S 
when  P=90  degrees.  Since  the  relation  My  =  Rcoty  Fx  /  4  holds  both  for  finite  S  and  in  the 
hyperthermal  limit,  plots  of  the  relative  error  in  Mv  are  identical  to  those  in  Figure  5.8. 

The  relative  error  in  the  hyperthermal  approximation  toAf.(fi))  is  shown  in  Figure  5.9  for  |3=0,  45,  & 
90  degrees.  In  contrast  with  My(Fx), the  relative  error  now  increases  with  increasing  angle  of  attack 
for  a  given  molecular  speed  ratio. 


Relative  Error  in  Fx 


amz  Relative  Error  in  M- 


Figure  5.9.  Cone  Hyperthermal  Spin  Moment  Error 
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Appendix:  Elements  of  the  Kinetic  Theory  of  Gases 

The  purpose  of  this  appendix  is  to  present  a  concise  account  of  those  aspects  of  gas  kinetic  theory 
that  are  utilized  in  the  report  to  calculate  the  aerodynamic  disturbances  on  a  body  immersed  in  a 
free-molecular  flow.  We  restrict  our  attention  to  ideal  gases  and  develop  the  simple  equation  of 
state  and  the  Maxwell-Boltzmannn  velocity  distribution  from  considerations  of  molecular 
collisions.  Further  details  can  be  found  in  [3]  &  [4]. 


The  Ideal  Gas 

It  proves  convenient  to  define  an  ideal  gas  by  imposing  the  following  assumptions  which  have 

been  shown  to  be  extremely  reasonable  in  many  areas  of  application. 

1)  The  gas  is  composed  of  identical  molecules  which  are  separated  by  distances  which  are  large 
compared  with  the  molecular  diameter.  Any  finite  volume  of  gas  contains  an  extremely  large 
number  of  molecules.  The  first  assumption  will  be  valid  for  gases  under  low  pressure.  To 
appreciate  the  significance  of  the  second  assumption,  we  recall  the  experimental  result  that  a 
1  Kgm-mole  of  gas  contains  approximately  6.02214  xlO26  molecules  (Avogadro’s  number). 
Since  at  standard  conditions.  1  Kgm-mole  of  gas  occupies  a  volume  of  22.4  m\  we  see  that 
there  are  approximately  3x  10'°  molecules  in  a  cubic  millimeter  at  standard  conditions. 

2)  In  the  absence  of  external  forces,  the  molecules  are  distributed  uniformly  throughout  their 
container.  Thus  if  their  are  a  total  of  N  molecules  in  a  container  of  volume  V.  then  the  average 
number  of  molecules  per  unit  volume,  n,  is  n  =  N/V .  It  follows  that  a  small  volume  element 
AV  contains  n  AV  molecules.  Clearly,  it  is  possible  that  by  taking  AV  sufficiently  small  no 
molecules  will  be  contained  in  the  volume  element.  However,  since  the  number  density  of 
molecules  is  so  large,  we  can  divide  the  container  into  extremely  small  elements  of  volume 
which  still  contain  a  large  number  of  molecules.  Even  a  cube  1/1000  mm  on  a  side  will 
contain  approximately  3xl07  molecules  (at  standard  conditions).  This  dimension  is 
infinitesimal  compared  with  other  physical  dimensions  of  interest.  We  therefore  are  justified 
in  applying  the  methods  of  differential  and  integral  calculus  to  these  “small”  volume  elements 
with  negligible  error.  Thus  the  number  of  molecules  in  a  differential  volume  element  dV  is 
given  ndV . 

3)  Forces  are  exerted  upon  a  molecule  only  when  it  collides  with  another  molecule  or  the  walls 
of  its  container;  these  collisions  are  regarded  as  perfectly  elastic.  The  walls  of  the  container 
are  assumed  to  be  smooth  hence  there  will  be  no  change  in  the  tangential  component  of 
molecular  velocity  upon  collision. 

4)  There  is  no  preferred  direction  regarding  the  molecular  velocity.  Imagine  a  vector  attached  to 
each  molecule  with  magnitude  and  direction  scaled  to  its  velocity.  Transfer  these  vectors  to  a 
common  origin  and  construct  a  sphere  of  arbitrary  radius  r  centered  at  the  origin.  Our 
assumption  implies  that  these  vectors  (prolonged  if  necessary)  will  intersect  the  sphere  at  a 
set  of  points  uniformly  distributed  over  its  surface.  The  number  of  these  points  will  be  N,  the 
number  of  molecules,  and  their  surface  density  is  given  by  N/^Anr2}.  Introducing  spherical 

coordinates  { r,0,(p }  at  the  center  of  the  sphere  where  0  is  the  polar  angle,  we  find  that  the 
number  of  molecules  per  unit  volume  with  velocities  having  directions  between  (6,6 +  d6) 
and  ( (p,(p  +  d(p )  is  given  by 


56 


Yl 

dnfkf  =  —  sin  8  ddd(p  (A.  1) 

5)  Clearly,  the  speed  of  any  given  molecule  changes  continually  due  to  collisions  with  other 
molecules  or  the  walls  of  the  container.  To  simplify  the  mathematics,  we  assume  that  the 
molecular  speed  can  assume  any  value  in  the  range  (0,°°)  and  that  the  number  of  molecules 
with  speeds  in  a  specified  range  is  constant  for  a  gas  in  thermal  equilibrium.  Rigorously 
speaking,  there  will  be  a  finite  upper  bound  on  the  molecular  speed.  Later  we  will  derive  a 
statistical  distribution  of  molecular  speeds  (based  on  these  assumption)  and  show  that  this 
distribution  function  tails  off  very  rapidly  predicting  a  very  small  fraction  of  molecules  with 
speeds  in  a  high  range.  This  provides  justification  for  removing  the  finite  upper  bound  on  the 
molecular  speed. 

Molecular  Surface  Collisions 

It  is  important  to  consider  the  number  of  molecular  collisions  per  unit  time  per  unit  area  of 
material  surface.  Figure  A1  shows  an  element  of  surface  area  dA  with  normal  n  along  with  the 
polar  and  azimuth  angles  6,(p.  Molecules  impinge  upon  the  surface  from  all  directions  with  all 
speeds.  Let  us  consider  only  those  molecules  traveling  with  specific  speed  v  and  specific  direction 
6,(p  which  we  term  a  “  Otpv  -molecule”.  More  precisely,  these  molecules  have  speeds  between  v 
and  v+Av,  and  directions  in  the  range  (6,0 +  dd)  and  (cp,tp  +  d(p).  In  a  time  interval  dt,  the 
molecules  which  impinge  on  the  surface  are  contained  in  a  cylinder  of  slant  height  vdt  and  base 
area  dA  with  corresponding  volume  vcosOdtdA.  If  we  denote  by  dnv  the  number  of  molecules 
per  unit  volume  with  speeds  between  v  and  v+dv,  then  dnjn  represents  the  fraction  of  molecules 
in  this  speed  range.  We  make  the  further  assumption  that  this  fraction  also  applies  to  any 
subgroup  of  molecules  e.g.  those  traveling  in  any  particular  direction.  Utilizing  eq.(A.l),  we  find 
that  the  number  of  9<pv  -molecules  per  unit  volume  is  given  by 

— sin  9  d6  dm  dnv 
4  n 

Thus  the  number  of  6(pv  -molecules  in  the  cylinder  i.e.  the  number  of  collisions  in  time  dt  is 

—  vsm  9  cosddO  dm  dnrlAdt  (A. 2) 

4  n 

and  the  number  impinging  per  unit  time  per  unit  area  is  given  by 

—  v  sin  9  cosOdO  dtp  dn 
An 

By  integrating  the  above  expression  over  \0<9  <n/2\  0<(p<2n },  we  obtain  the  number  of 
molecules  with  speed  v  colliding  with  the  surface  from  all  directions  per  unit  time  per  unit  area 


In  accordance  with  assumption  (5)  above,  the  total  number  of  molecules  impinging  upon  the 
surface  per  unit  time  per  unit  area  is  given  by 
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Figure  A 1  Molecule-Surface  Collision 


jfvr/n,.  (A.3) 

4Jo 

We  next  relate  the  above  integral  to  the  average  molecular  speed.  First  consider  the  discrete  case 
of  N  molecules  for  which  N]  have  the  same  speed  v,,  N2  have  the  same  speed  v2,  etc.  The  average 
speed  v  is  then  given  by 

-=Nivl+N2v2+... 

N 

Dividing  numerator  and  denominator  by  the  volume  V  of  the  container,  we  can  write 

_  n.v,  +/j,v,  +... 

V  =  — - — - 

n 

where  ni  is  the  number  of  molecules  per  unit  volume  with  speed  v,.  Transitioning  to  the  case  of  a 
continuous  velocity  distribution,  we  have 

v  =  —  [  vdnv 
n  0 

Utilizing  this  result  in  eq.(A.3)  we  arrive  at  the  desired  expression  for  the  total  number  of 
molecules  impinging  upon  the  surface  per  unit  time  per  unit  area 

-«v  (A.4) 

4 

We  now  proceed  to  calculate  the  force  imparted  by  the  impinging  molecules  upon  the  surface. 
With  reference  to  Fig.  Al,  consider  a  single  Qcpv  molecule  of  mass  m  before  and  after  a  surface 
collision.  Regarding  the  material  surface  to  be  of  infinite  mass  relative  to  the  molecules,  and 
recalling  assumption  (3)  above,  we  find  the  change  in  momentum  per  molecular  collision  to  be 
2mvcos0n.  Multiplying  by  (A. 2).  we  obtain  the  change  in  momentum  dp  for  all  8(p\>  molecules 
colliding  with  the  surface  in  the  time  interval  dt. 

dp  =  —  niv2  sin0cos20  ddd(pdnvdAdtn 
2  n 
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Integrating  this  result  over  the  region  {0<9<7i/2,  Q<(p<2n,  0<v<°°}  and  dividing  by  dt 
gives  the  average  force  exerted  by  the  surface  on  all  the  incident  molecules.  By  Newton’s  third 
law,  the  force  exerted  by  the  molecules  on  the  surface  are  oppositely  directed  with  an  associated 
pressure  p  given  by  the  expression 


Following  the  same  argument  as  given  above  for  the  average  speed,  we  can  express  this  integral 
in  terms  of  v2  -  the  average  value  of  the  square  of  the  speed. 


where  p  denotes  the  gas  density. 


1  —  1  - 
p  =  —nmv  =—pv 
3  3 


(A.5) 


The  above  relation  can  be  used  to  arrive  at  a  molecular  interpretation  of  absolute  temperature. 
The  equation  of  state  of  an  ideal  gas  can  be  expressed  in  the  form 


pV^nKT 

where  Ti  is  the  number  of  moles,  7?  is  the  Universal  Gas  Constant  (8314.5  joules /Kg -mole  °K), 
and  T  is  the  absolute  temperature.  But  the  number  of  moles  is  equal  to  the  total  number  of 
molecules,  N,  divided  by  the  number  of  molecules  per  mole,  N0-Avogadro’s  number 
(6.02214 xlO26  molecules /Kg -mole).  The  universal  constant  <R/N0  =  k  is  Boltzmann’s 

constant  (1.38065  xlO-23  joules/molecule  °K),  in  terms  of  which  the  ideal  gas  law  assumes  the 
form 


pV  =  NkT  (A.6) 

Since  n  is  the  number  of  molecules  per  unit  volume,  N/V,  eq.(A.5)  can  be  expressed  in  the  form 


pV  -  —  Nmv 
3 


If  this  is  to  agree  with  eq.(A.6),  we  must  conclude  that 


(A.7) 


which  says  that  the  mean  translational  kinetic  energy  of  a  gas  molecule  is  proportional  to  the 
absolute  temperature. 
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Distribution  of  Molecular  Velocities 


In  this  section  we  derive  the  celebrated  Maxwell-Boltzmann  velocity  and  speed  distributions  for 
an  ideal  gas.  A  rigorous  derivation  requires  methods  from  statistical  mechanics  which  will  not  be 
attempted  here.  We  present  essentially  the  argument  originally  given  by  Maxwell  which  although 
heuristic,  is  plausible  from  a  physical  point  of  view.  In  the  previous  section  we  arrived  at 
expressions  for  the  molecular  flux  and  pressure  in  terms  of  the  mean  and  rms  molecular  speeds 
respectively.  By  imposing  additional  assumptions,  we  can  calculate  the  actual  velocity 
distribution  which  will  ultimately  lead  to  relations  between  these  average  velocities  and  the  gas 
temperature. 


Consider  a  rectangular  coordinate  system  whose  axes  {y,.,vrv,}  correspond  to  the  three 
components  of  molecular  velocity.  Each  molecule  in  physical  space  has  a  corresponding  point  in 
this  “velocity  space”.  Since  the  number  of  molecules  if  finite,  the  distribution  of  these  points  in 
the  velocity  space  will  not  be  continuous.  However,  in  view  of  assumptions  (1)  &  (2)  above,  we 
will  approximate  this  distribution  by  a  continuous  function  with  the  understanding  that  any 
designated  “volume  element”  in  the  velocity  space  must  be  chosen  large  enough  so  that  it 
contains  a  large  number  of  representative  points.  We  start  by  posing  the  question:  how  many 
molecules  (out  of  the  total  number  N)  have  x  components  of  velocity  between  y,  and  vx+dvx  ?  This 
is  equivalent  to  finding  the  number  of  points  in  velocity  space  contained  in  the  thin  slice  parallel 
to  the  v,-v.  plane  between  y,  and  vx+dvr  If  we  denote  this  number  by  dNv ,  then  the  fraction  of  the 

total  number  of  molecules  with  x  components  of  velocity  between  y,  and  vx+dvx  (regardless  of 
there  y  or  z  velocity  components)  is  dNv  /N .  This  fraction  is  certainly  a  function  of  v,.  It  seems 

plausible  that  it  should  also  be  proportional  to  the  thickness  of  the  slice.  Hence  we  write 


dN,. 

~N~ 


f(vMv, 


(A- 8) 


where  the  function /is  to  be  determined.  By  assumption  (4)  the  velocity  distribution  is  isotropic, 
so  we  must  have  for  the  same  function / 


7T  =  /(■’,)*. 

IN,.  .  . 

r=/(v')‘,v= 


(A.9) 


where  dNr  and  dN,..  are  defined  analogously  to  dNv  We  now  impose  the  additional  assumption 
that  these  fractions  are  also  applicable  to  any  subgroup  of  molecules  e.g.  for  the  subgroup  of 
molecules  with  x  components  of  velocity  between  vx  and  vx+dvx  f(v^dvy  of  them  will  also  have  a 
y  component  of  velocity  between  vv  and  vv  +  dv,..  Thus  the  number  of  molecules  simultaneously 
having  x  components  of  velocity  between  y,  and  vt  +  dvy,  y  components  between  v  and  vy  +  dvy 
and  z  components  between  v.  and  v.  +  dv.  is  given  by 

=  Nf{v>  )/(  Vr  )/(  v-  )dv,  dvy  dv: 

Geometrically,  this  gives  the  number  of  representative  points  in  the  volume  element  dvydvydv.  at 
the  location  (vt,vv,vj  in  velocity  space.  The  number  of  representative  points  per  unit  volume 
then  defines  a  “density”  p 
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P=W/(vJ/(vv)/(v_) 


The  change  in  p  for  small  increments  in  v^v  and  v.  is  given  by 

dp  =  ^-dv  +^~  dv  +  dv. 
dv,.  dv„  dv_ 


(A.  10) 


Since  the  velocity  distribution  is  isotropic,  there  will  be  no  change  in  p  if  we  restrict  ourselves  to 
those  points  that  are  at  a  fixed  distance  from  the  origin  in  velocity  space  i.e.  dp-  0  under  the 
constraint  v2  +  v2  +  v2  =  const.  Using  eq.(A.lO),  this  results  can  be  expressed  in  the  form 


lMdv  +ZfcL, 


dv„  +  - 


/M  '  f(v,)  ’  fM  ' 

for  velocity  increments  satisfying 


dv ,  =  0 


vYdvx  +  vydvr  +  v.dv,  —  0 


(A.  11) 


(A.  12) 


In  eq.(A.ll)  only  two  of  the  differentials  are  independent.  Utilizing  eq.(A.12)  we  can  recast 
eq.(A.ll)  in  the  form 


f\vs) 

JM  fMvz 

from  which  it  follows  that 


dvx  + 


f(v>) 


fW 

/(v-)  ' 


dv„ 


0 


/'(-■)  /'(”■)  AP 

v,/(v.)  v,f(v,)  v./(v.) 


The  first  ratio  is  a  function  only  of  v(,  while  the  second  and  third  are  respectively  functions  only 
of  vv  and  v_.  It  follows  that  each  is  equal  to  a  constant 


f'M  _  i 

v,/W“ 


etc. 


Integrating  this  equation  and  observing  that  f(vx)  — >  0  as  vv  — >  we  obtain 


(/32  -  A/2) 


/(u)  =  aexp(-/3V) 


(A.  13) 


Thus  we  have  obtained  the  desired  form  of  the  velocity  distribution  function/.  In  the  next  section 
we  will  obtain  the  physical  interpretation  of  the  parameters  a  and  (5 .  We  can  now  express  the 
number  of  molecules  with  x  components  of  velocity  between  yv  and  vx+dvx  in  the  form  (see 
eq.(A.8)) 


dNv  -  Na  exp(-/32v2 ) dvx  (-°o  <  vr  <  °o) 

(with  similar  expressions  for  dNv  and  dNv )  and  the  number  of  points  per  unit  volume  (A.  10)  as 
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p  =  /Va’exp^jSV) 


(A. 14) 


where  v  is  the  molecular  speed. 

Many  times  we  are  interested  in  the  number  of  molecules  with  speeds  in  a  certain  range  without 
regard  to  direction.  The  speed  range  of  interest  may  be  of  finite  or  infinite  length.  Let  us  first  find 
the  number  of  molecules  with  speeds  in  the  infinitesimal  range  between  v  and  v  +  dv.  Denoting 
this  number  by  dNr,  we  simply  need  to  find  the  number  of  representative  points  inside  the  thin 
spherical  shell  of  radius  v  and  thickness  dv  in  velocity  space.  Since  p  is  uniform  in  this  region, 
we  simply  multiply  the  density  given  in  eq.(A.14)  by  the  volume  4 nv2dv  of  the  shell 

dNr  =  4 nNa*v2  exp (A.  15) 

Note  that  the  velocity  distribution  (A.  14)  assumes  its  maximum  value  at  v=0  and  decays 
exponentially  to  zero,  while  the  speed  distribution  (A. 15)  is  zero  at  v=0.  increases  to  its  maximum 
value  and  then  decays  exponentially  to  zero.  The  number  of  molecules  with  velocities  in  a  desired 
range  or  having  speeds  in  a  specified  interval  can  be  calculated  by  integrating  the  appropriate 
distribution  over  the  appropriate  domain.  Thus  the  number  of  molecules  with  speeds  between  v, 
and  v,  (0<  v,  <  i'2)  is  obtained  by  integrating  (A.  15)  over  the  interval  v,<v<v2. 


Evaluation  of  the  Distribution  Parameters 

We  now  show  how  the  two  velocity  distribution  parameters  a  and  P  can  be  related  to  the  absolute 
temperature  of  the  gas.  Firstly,  the  integral  of  (A.  15)  over  all  speeds  from  zero  to  infinity  must 
yield  TV-the  total  number  of  molecules.  Performing  the  integration,  we  find  that 

a  =  P/Jrt 

The  speed  distribution  (A.  15)  can  then  be  expressed  in  terms  of  the  single  parameter  P 

dNy  =  pV  exp(-j 32v2}dv 

We  now  relate  the  parameter  P  to  the  most  probable,  mean,  and  rms  speeds.  We  define  the  most 
probable  speed  v((  to  the  value  of  v  for  which  dNJdv  assumes  its  maximum  value.  It  is  easy  to 
show  that 

v„  =  P~' .  (A.  16) 


The  mean  speed  v  is  given  by 


The  rms  speed  is  defined  by 


(A.  17) 


Combining  this  last  relation  with  eq.(A.7),  we  deduce  that 

P  =  jm/2kT 


(A.  18) 


allowing  us  to  express  all  the  distribution  functions  in  terms  of  absolute  temperature  and 
molecular  mass. 


